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AN EVALUATION OF RANDOM ANALYSIS METHODS 


FOR THE DETERMINATION OF PANEL DAMPING 

By Waman V. Bhat and John F. Wilby 
The Boeing Company, Seattle, Washington 


1.0 SUMMARY 


An analytical study is made of steady-state and non-steady-state methods for the mea- 
surement of panel damping. Particular emphasis is placed on the use of random process tech- 
niques in conjunction with digital data reduction methods. The steady-state methods considered 
use the response power spectral density, response autocorrelation, excitation-response cross- 
power spectral density, or single-sided Fourier transform (SSFT) of the response autocorrela- 
tion function. Non-steady-state methods are associated mainly with the use of rapid frequency 
sweep excitation. Problems associated with the practical application of each method are 
evaluated with specific reference to the case of a panel exposed to a turbulent airflow, and 
two methods, the power spectral density and the single-sided Fourier transform methods, are 
selected as being the most suitable. These two methods are demonstrated experimentally, and 
it is shown that the power spectral density method is satisfactory under most conditions, pro- 
vided that appropriate corrections are applied to account for filter bandwidth and background 
noise errors. Thus, the response power spectral density method is recommended for the mea- 
surement of the damping of panels exposed to a moving airflow. 



2.0 INTRODUCTION 


The amplitude of vibration of a structural system exposed to an exciting force is con- 
trolled, to a large extent, by the total damping present in the system. When predicting the 
response of a structure to a given excitation it is necessary to know the damping to a suffi- 
cient degree of accuracy. Equally important, when measuring the vibration of a model or 
full-scale structure under certain excitation conditions, measurements of the panel damping 
are required if the experimental data are to be extrapolated to other conditions with reason- 
able reliability. 

Techniques for the measurement of the damping of a structure have been developed 
over the years, but, until recently, the methods have been limited to the use of steady-state, 
discrete frequency excitation (e.g., refs. 1 through 4). For a structural system whose damping 
does not vary with operating condition or environment, and where measurements can be 
made under specially controlled test conditions, the discrete frequency methods are usually 
adequate. However, when the damping is dependent on operating condition, for example 
where flutter or acoustic damping is important, alternative techniques for damping measure- 
ment have to be used. In response to this need, techniques using steady-state random excita- 
tion and transient excitation have been developed in recent years. These alternative methods 
can be used when the structure is exposed to random excitation or is exposed to the operat- 
ing environment for only a short time or when noise interference occurs in the measured 
excitation or response. 

The case of a structure exposed to a moving airstream is an important example of a test 
situation where the steady-state, discrete frequency methods are unreliable. It was for this 
reason that the present study of damping measurement methods was undertaken. Results of 
the study, which considers steady-state random excitation and transient excitation methods, 
are contained in this report. Essentially, the discussion can be separated into two parts, the 
first of which is an analytical survey of available methods. The second part describes an 
experimental demonstration of two techniques that, on the basis of the analytical study, are 
judged to be the most suitable for the application under consideration, namely, the measure- 
ment of panel damping in the presence of a moving airstream. Throughout the discussion, 
emphasis is placed on the accuracy and reliability of each method in practical situations. 
Digital techniques for data reduction will be required for several of the damping measurement 
methods. Therefore, special attention is given to problems associated with these digital 
techniques. 

As a final result of the study, techniques will be recommended for the measurement of 
panel damping in the presence of an airflow. In principle, the methods will be applicable to 
structures of greater complexity than a simple panel in a wind tunnel wall. However, it may 
be necessary to assume that there is some knowledge of the natural frequencies and mode 
shapes of the system so that the measuring tranducers can be placed at optimum locations. 

The program was carried out at The Boeing Company, Seattle, Washington and was 
monitored by Dr. Wei J. Chyu, NASA Ames Research Center. 
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area of structure, for panel A = L j L 3 
see equation (29) 

defined by equations (43) and (83) 
defined by equation (84) 

Viscous damping rate 

real part of cross-power spectral density function Sp (x x ”, 00) 

real part of cross-power spectral density function Sp] (x ', x ”, co) 

generalized damping coefficient in mode 

equivalent viscous damping coefficient in terms of the hysteretic 
damping, (k^Vcj 

flexural rigidity of a plate 

weighting function, lag window 

see equation (39) 

single-sided Fourier transform of the response autocorrelation 
function, R w (r), defined by equation(71) 

frequency 

center frequency of analyzing filter 

natural frequency of mode (a = m,n) 

upper frequency of interest 

ideal half-power bandwidth of a resonance peak 

filter bandwidth at the half-power point 

measured bandwidth of a resonance peak 

frequency interval at which frequency-dependent functions are 
calculated, 1/AT' 

modulus of vector H 

complex frequency response function for mode of order a, 
defined by equations (4) and (5) 




complex conjugate of H a (cj) 

H a (to, t) 

time-varying frequency response function 

lH a (co)i 2 

time-averaged square of frequency response function 

h a (t), h a (r) 

impulse response function for mode of order a, defined by 
equation (5) 


defined by equation (18) 

i 

nAT 

j 

positive integer 


displacement joint acceptance for mode a 

V M > 

displacement cross acceptance for modes a , (3 

K, K', K", K a , 

miscellaneous constants or functions, defined in text 

k 

time constant in co '(t) 

L 

number of degrees of freedom in section 4.1.7 

L 

number of data points used for calculating finite discrete 
transforms, 2J 

L 1 

panel length in xj direction 

L 3 

panel length in x^ direction 

L« 

generalized external force in mode of order a 

M 

mass per unit area 

M a 

generalized mass in mode of order a 

m 

mode order in x j direction denoting the number of modal 
half-wavelengths 

N 

number of measurement locations in section 4.1.7 

N 

number of data points 

n 

mode order in X3 direction, denoting the number of modal 
half-wavelengths 

n(t) 

noise in response (see fig. 1 ) 

n 

ratio of ideal bandwidth of resonance peak to filter bandwidth. 
Af a /Af F 



P(x,co) 

P*(x,CO) 

Pj(to), P2(co) 
P(C4 

P = p(x» t) 
Pl(x'»t)> P2(i"> t ) 
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Q p (x', x", co) 

Qpl 
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q a (t) 
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R w (x» T ) 
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R W ( T )> R n( T ) 
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S(w) 

Sp(oj) 

S p] (oj), S p2 (oj) 


Fourier transform of p(x,t), defined by equation (10) 
complex conjugate of P(x,a>) 

Fourier transform of pj(t) and p 2 (t), respectively 

Fourier transform of p(t) 

pressure at point x and time t 

excitation and noise in excitation (see fig. 1) 

pressure at time n • At 

complex part of Sp (x', x", co) 

complex part of Sp](x', x", c_>) 

Fourier transform of q a (t) 
generalized coordinate of mode of order a 
autocorrelation function 
excitation pressure cross-correlation function 
displacement autocorrelation function 
narrow band autocorrelation function 
defined by equation (49) 
defined by equation (5 1 ) 
cross-correlation function, equation (62) 
autocorrelation at time delay n • At 
mode order in Xj direction 

defined by equation (39) 

white noise power spectral density 

power spectral density function 

power spectral density function of excitation p(x, t) 

power spectral density functions of excitations pj(x, t) and 
p 2 (x, t), respectively 



S w (x, w) 

S w (x, w) a 
S w p (x, x'> w) 

S w p (x, x',w) a 
S k (k,Af F ) 
s 
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ds/df 
t, T 
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X a 

X = (Xj, X 3 ) 

a = (m,n) 

0 — (r,s) 

^pn 

^pn 

r mn^) 
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power spectral density function of the displacement w(x, t) 

displacement power spectral density function in mode a 

cross-power spectral density function relating the excitation 
p(x', t) and response w(x, t), equations (54) and (55) 

contribution of a mode to S w p(x, x', a>) 

power spectral density at a frequency of k • Afp 

digitizing or sampling rate, 1/Ar 

mode order in x 3 direction 

rate of change of arc length with frequency 

time parameters 

total sample length 

length of sample subinterval, N • Ar 

length of subinterval after addition of zeros 

Fourier transform of L a (t) 

number of subspectra, T/AT 

Fourier transform of the displacement w(x, t) 

displacement at point x, and time t 

see equation (35) 

distance vector in (xj ,x 3 ) plane 

two-dimensional mode order 

two-dimensional mode order 

defined by equation (34) 

complex conjugate of (3p n 

defined by equation (3 1 ) 

nondimensional cross-power spectral density function relating 
the excitation and response, defined by equation (56) 

1 for p = q, 0 for p# q 
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delta function 1 for x' = x 0 , 0 for x'#x Q 

delta function, 1 for r = 0, 0 for r ^=0 

normalized standard error, equation (92) 

error in the frequency response function due to the drift in 
natural frequency 

viscous damping factor in mode of order a 
equivalent viscous damping factor 
see equation (25) 

sum of response power spectral densities at all measuring locations 
due to mode of order p, defined by equation (36) 

see equation (39) 

stiffness 

stiffness in the mode of order a 
unitary matrix, defined by equation (32) 
complex conjugate of Apj 
panel mass per unit area 

hysteretic damping factor for mode of order a 

radius of curvature of the cross-power spectral density curve 

ratio of background noise power spectral density to the true 
response power spectral density at resonance 

time delay 

maximum time delay in autocorrelation function 
sampling interval, 1/s 

spectral window, filter frequency response function 

time-varying phase angle 

defined by equation (69) 

panel mode shapes of order a. and |3 
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Aco 
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copplex conjugate of \p a and respectively,. = *P a and 

H S H 

angular frequency 

angular frequency of mode of order a 

, ? J/2 

natural frequency for damped free vibrations, co a (1 - 

drift in the natural frequency 

maximum value of the drift in natural frequency 

time-independent part of resonance frequency of a system with 
a time-varying frequency response function 

angular frequency bandwidth 

angular frequency half-power bandwidth of a resonance peak 
for mode of order a 



3.0 DAMPING OF STRUCTURES 


3.1 GENERAL DISCUSSION 

Two important sources of damping in typical aerospace structures are structural hysteresis 
and acoustic radiation. For structures with no special damping treatments, hysteretic damping 
is mainly due to the joints in the structure, although there is some contribution from the 
damping of the material itself. The hysteretic damping force is proportional to displacement 
and is in counterphase with velocity. Acoustic radiation damping arises from the dissipation 
of energy in the surrounding air and is viscous controlled. Viscous damping forces are propor- 
tional to, and in counterphase with, the structural velocity. 

From available literature it is possible to determine typical values of the damping ratios 
of structures ranging from single rectangular plates to stiffened airplane structures. Maestrello 
(ref. 5) and Wilby (ref. 6) show damping ratios in the range 0.001 < f < 0.008 for single 
rectangular panels with clamped boundaries. In some cases (ref. 6), the damping is seen to be 
dependent on the environment, presumably when the acoustic radiation damping forms a 
significant fraction of the total damping. Laboratory measurements on stiffened panels by 
Mead (refs. 7 and 8) give damping ratios in the range 0.002 to 0.013. Higher damping ratios, 
0.014 to 0.020, have been measured on aircraft structures by Clarkson and Ford (ref. 9), 
although Mead (ref. 7) has placed some doubt on the accuracy of these values because the 
method of data anlaysis could cause overestimation of the damping. 

The damping of structures can be increased greatly by the addition of special treatments 
such as damping tape. Mead (ref. 8) lists damping ratios up to 0.07 for a rectangular aluminum 
plate with one layer of damping tape, and Maestrello (ref. 5) measures damping ratios as high 
as 0. 1 for a plate with two layers of tape. These results show a more than ten-fold increase in 
the damping with respect to the untreated case. The damping introduced by the treatment, 
using either unconstrained or constrained layers, is hysteretic in nature (ref. 8), with the 
result that the damping of the treated structure will be mainly hysteretic. 

The above data show that, for untreated structures of the type likely to be found in 
aerospace construction, the damping will be low. Thus, in the mathematical analysis, the 
damping can be represented for simplicity as either solely hysteretic or solely viscous. As 
will be seen later, this is particularly useful for random vibration where the concept of 
hysteretic damping is not strictly valid. When damping treatment is added and the damping 
is mainly hysteretic, the assumption that the damping is solely viscous will still be made, 
although it will introduce some error. 

Several methods for measuring the damping of structures are available— discrete frequency 
or random, steady-state or transient— and the choice of method for a particular test will 
depend on the experimental conditions present. Discrete frequency methods are suitable for 
conditions where there is no noise present in either the excitation or response. These methods 
have been discussed extensively by several authors (see refs. 1, 2, 4, and 10) and will not be 
included in this report. In the present case of interest, the steady-state, discrete frequency 
method is unsatisfactory because of the high level of “noise” introduced by the pressure 
fluctuations in the airflow, and the method has been excluded from consideration. Thus, the 



analytical study in sections 4.0 through 6.0 will consider only steady-state random methods 
and transient methods. In some cases the turbulent airflow will provide the necessary excita- 
tion, but in other cases a second excitation will be introduced as a reference and the turbu- 
lent airflow regarded as noise in the system. The methods described will apply to any repre- 
sentation of the panel vibration (displacement, velocity, acceleration, etc.), but since the 
method chosen will be used with displacement probe instrumentation, the structural vibra- 
tion will be described, without loss of generality, in terms of displacement. 

Several of the measurement methods discussed will be suitable for either analog or 
digital data analysis techniques, but in other cases digital analysis procedures will be neces- 
sary. Thus, emphasis will be placed on the use of digital techniques for all methods, particu- 
larly when these techniques involve relatively new problem areas. 


3.2 PANEL VIBRATION 

At several stages in the analysis of the damping methods, reference will be made to the 
equation of motion for rectangular panels. Thus it is appropriate to outline the normal 
mode analysis of a rectangular plate exposed to random excitation, the resulting displace- 
ment power spectral density and autocorrelation equations being extensions of the corre- 
sponding equations for single-degree-of-freedom systems. Although the analysis is performed 
with specific reference to rectangular plates, the general conclusions will have wider application. 

Assume that the displacement w(x,t) of the vibrating structure obeys the thin-plate 
equation for small deflections 


^ + C j7 + DV 4 w = p(x,t), (1) 

dt 2 dt 

where M is the mass per unit area, C is the viscous damping, D is the flexural rigidity, and 
p(x,t) is the exciting force. D can be complex if hysteretic damping is included in the stiffness 
term. Assume further that the solution to equation (1) can be expressed in terms of a linear 
sum of normal modes, i.e., 


W(x,t) = £ q (t)t|/ (x), (2) 

a 

where q^(t) is the generalized coordinate for the ath mode and i// a (x) is the mode shape 
function. Then the equation of motion for the ath mode is given by 

where M a , C a , /c a (l + iy a ), and L a are, respectively, the generalized mass, viscous damping, 
complex stiffness, and external force in the oThmode. The parameter v a represents the 
hysteretic damping of the system. If the panel surface density is denoted by /x(x), then 

M a = / / i Q0'/' a 2 (x)dx, 

A 



where the notation J ^ dx denotes the double integral Jq Jq dx3dx j over the panel 
areaA = LjL3« 

For purely viscous damping the frequency response function for the mode is 

H a (oS) = ~ \ ’ (4) 

M a [(co a Z -co 2 ) + 2i^0J a 0J] 

where co a = y/K a /M a is the undamped natural frequency and = C a /2VA^K^is the damping 
ratio. The frequency response function H a (o;) and the impulse response function h a (t) are 
related by Fourier transformation, such that 


h <*(t) = 27jf*H a Me lwt doj 

and (5) 

H a M = j[~ h a (t)e“ iwt dt, 

it being noted that h a (t) = 0 for t < 0. It can be shown that h a (t) is real, as is expected from 
physical reasoning. 

For hysteretic damping alone, 


H a M - 


1 t 

M a [ (co a 2 - co 2 ) + i v a ^a ] 


( 6 ) 


but the Fourier transform of this function is complex and does not represent a physically 
realizable response. This inconsistency illustrates the difficulties encountered when the con- 
cept of hysteretic damping is used in the analysis of vibrations which are not of the deter- 
ministic sinusoidal type. It has been suggested (ref. 8) that, to overcome the problem, the 
hysteretic damping term be replaced by an equivalent viscous term C ae = (K a ^ a )/oj. The 
equation of motion would then be 




• (7) 


Then the analysis can be performed as for viscous damping and the hysteretic damping rein- 
stated in the final result (replacing the term tof ae /co a with v a /2). This argument provides 
further justification for considering only viscous damping in the subsequent analyses. 

If the equation of motion for viscous damping, similar to equation (7), is Fourier trans- 
formed with respect to time, the equation reduces to 


Q a (co) = H a (cu)U a M 


( 8 ) 



where Q^Cco) is the Fourier transform (or complex Fourier amplitude) of q a (t), i.e., 

Q a (w)= ijC q « (t)e_ia;tdt ’ 

and U a (co) is the Fourier transform of L a (t). 

Also, 


U a (w) = 2~ J [ oo e_1 ° Jt P(x,t)i// a (x)dxdt 

= J^t// a (x)P(x,uj)dx 

where 

P(x,co) = jf ^p(x,t)e~ 10jt dt. 


(9) 


( 10 ) 


Thus, from equation (2) the Fourier transform of the displacement is 

W(x,<4 = £Q>j)«//Jx), (11) 

a u w 

and the displacement power spectral density function at x = (x is 

= W(x,co)W*(x,co), (12) 

where the asterisk denotes the complex conjugate. The parameter x in S w (x,gj) is included to 
indicate that the displacement power spectral density function for a finite structure is 
inhomogeneous. 

Substituting equations (8), (9), and (1 1) in equation (12), 



where Sp(x',x".co) is the excitation pressure cross-power spectral density function, 

S p (x',x",co) = P(x',Go)PV,co). 


Using the relationships 


R(t) = r°°S(oj)e i6jr doJ 

-oo 


S(^) = ~f\(T) C -^ T dT , 


and 



the displacement autocorrelation function at x is 


R w (x,r) = S A 4^ a (x')^(x" )X*H a *(w)H^(w) S p (x',x‘ w J w)e iwT dwdx’dx" 

“ 5 “ _ (14) 

= ¥|'*a ( - ) V 5) ii'''“ <!: ’ > ^ l!£ ’ )j C-C h “ <Tl)h ' S< ’’ 2) 

Rp (x',x", r + Tj - T 2 )dTjdT 2 dx'dx", 

where R p (x', x' , r) is the excitation pressure cross-correlation function. 

For a single-degree of freedom system such as a dashpot, the above equations (13) and 
(14) reduce to the much simpler forms 

S w (co) = |H(6j)| 2 S p (co) (15) 

and 

R w (r) = /°° f°° h(Tj )h(r 2 ) R D (r + tj - T 2 )drj dT 2 , ( 1 6) 

•7— OO*'-— OO * ~ 

respectively. 


3.3 FREQUENCY RESPONSE FUNCTIONS 

In some of the methods discussed in the following sections, the procedure reduces to the 
plotting of the frequency response function H a (co) in the complex plane. The damping ratio 
is estimated by fitting a circle to the data points in the neighborhood of the natural frequency, 
as a representation of the single-degree-of-freedom system. This technique is used to minimize 
effects which cause the measured function to deviate from the true frequency response func- 
tion for a single mode. 

Kennedy and Pancu (ref. 4) have shown that for a single-degree-of-freedom system with 
hysteretic damping, the plot of H~*(oj) is a straight line in the complex Argand plane and 
the plot of H a (co) (equation (6)) is a circle. However, for viscous damping the locus of H^(oo) 
is a parabola and can be considered as an approximation to a straight line only in the neighbor- 
hood of the natural frequency when the damping is small. Thus, the fitting of a circle to the 
experimental data for H a (co) is, in itself, an approximation when the damping is considered 
to be viscous. 



4.0 STEADY-STATE METHODS 


In this section, four methods using steady-state random excitation will be discussed. 

The methods use: 

a) The response power spectral density function 

b) The response autocorrelation function 

c) The excitation-response cross-power spectral density function 

d) The single-sided Fourier transform of the response autocorrelation function 

In three of the methods, a), b), and d), the pressure fluctuations in the turbulent airflow 
provide the random excitation forces. Method c) requires the introduction of a second excita- 
tion that is used as reference for the cross-power spectral density function, since it is not 
practical to use the airflow pressure fluctuations as a datum. 


4.1 RESPONSE POWER SPECTRAL DENSITY 


4,1.1 Introduction 

The response power spectral density method is the stochastic equivalent of the response 
amplitude method for discrete frequency excitation. The response power spectrum will 
exhibit a series of peaks associated with the natural frequencies of the structure, and the damp- 
ing, in common with the discrete frequency excitation case, can be estimated from measure- 
ments of the bandwidth of the resonant peaks. The response power spectral density method 
is conceptually simple since the damping factor is related to the bandwidth of the resonance 
peak. In principle, the method has the advantage that there are no requirements additional to 
those normally associated with the measurement of structural response (fig. la). However, 
for practical multimodal systems, the method suffers from a number of limitations that can 
be separated into two broad categories: first, the basic limitations inherent in the method 
and, second, the limitations imposed by the data analysis techniques and other peculiarities 
of the experimental procedure. The basic limitations include errors introduced by the close- 
ness of natural frequencies and the effect of statistical coupling. Examples of the second 
category will include errors due to drift of the natural frequencies, the finite bandwidth of 
the filter system, contributions from the background noise, and the statistical scatter due 
to the finite data sample length. These limitations are often interdependent and occasionally 
impose somewhat severe requirements on the experiment, e.g., if a very small filter bandwidth 
is required due to very light damping, then very long data sample lengths will be required to 
ensure acceptable statistical scatter (or to ensure high statistical reliability). 



Errors encountered in the practical use of the response power spectral density method 
are investigated analytically. The investigations are carried out with two questions in mind: 

a) How can the error under investigation be minimized? 

b) How can corrections be made to data acquired from a particular experiment? 

Some of the errors encountered are not peculiar to the power spectral density method 
since they, or their equivalents, arise in other methods. 


4.1.2 Basic Analysis 

The displacement power spectral density function for a rectangular panel exposed to 
random excitation is given by equation (13). It can now be written in the form 


S w (x,a>) =]£ £ x)0^(x)H tx * (co>H^ (co> ^(<^>» ( 1 7) 

where 

v w) = f A v- W)s P ( -'’- w,w)d -' d - w ( 1 8) 

and S w (x,w) is evaluated at point x = (xj.xg). Joint acceptances J aa (to) and cross acceptances 
J a p(oS) can be defined as 


W“) = 


4 

A 2 S p (co) 


I aa (w) 


and 




The joint acceptance term can be simplified slightly because if 
S p (x',x''^) = C p (x',x'',OJ) - iQ p (x',x",<4 

it can be shown (ref. 6) that 


J^J a ^'%(x")Q p ( x',x",w)dx' d x" = 0. 

It should be noted that the excitation p(x,t) will affect the response power spectral density 
S w (x,co) only through When the excitation is homogeneously distributed over the 



panel, the integrand in equation (18) is a function of separation distance only. For an 
inhomogeneous excitation, the integrand, and hence will depend on the location of 

the excitation on the panel. 

For a single mode, the displacement power spectral density is given by 
S w (x,w) = H a (co)H a *M^ a 2 (x)J aa (co)S p (cc) 


= |H a M| 2 ^ 2 (x)J aa (w)S p M. 

Thus, the displacement power spectral density depends on the joint acceptance J aa (co), loca- 
tion x (via mode shape ip a (x)), and frequency response function H a (oo). Jaa( w ) is weakly 
dependent on cu (see ref. 6) and can be considered to be constant in the vicinity of the reson- 
ance peak. The mode shape \fr a is essentially a scaling factor giving maximum deflection at 
an antinode and zero deflection at a node. Therefore, for a single mode, the displacement 
power spectral density is proportional to the square of the frequency response function, or 

S w (x,co) = K|H a (co)| 2 S p (w), 

where K is a factor that depends on the location of the measuring probe and the joint accept- 
ance. If the excitation pressure power spectral density is approximately constant in the fre- 
quency range of interest, the frequency response function and the damping of the system can 
be estimated from the displacement power spectral density function. The viscous damping 
ratio is given by 


_ AcUq, _ Afpt 

" 2to a “ 2f a ’ 


where Af a is the half-power bandwidth of the resonance peak of mode a, occurring at a 
frequency of f a . 

For a single-point excitation, with power spectral density S p (<j) at location x 0 , the 
displacement power spectral density function is 


S w (x,u>) 


S^ a 2 (x)|H a M| 2 J^^)// a (x0ij/ a (x'')S p M5(x' -x 0 )5(x w -x Q )dx'dx" 


E^ 2 (x)^ 2 (x 0 )|H a (co)| 2 S p (w). 


(19) 
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4. 1 .3 Statistical Coupling 


Equation ( 1 7) for the displacement power spectral density function can be written in 
the form 


S w (x,<4 =E'/' a 2 (x)|H a (co)| 2 l aa (w) 

oc oc p 

(«* 0 ) 


( 20 ) 


The cross terms, for which a & P, are due to statistical coupling between the modes of 
the panel, and these terms will cause the response peaks to deviate from the |H a (co)| 
shape. Thus, the statistical coupling or cross terms will influence the damping measurements. 


Three factors will affect the contributions from statistical coupling between modes oc 
and p: cross acceptances J a /j(co), frequency response functions H a (ou), and mode 

shapes t// a (x), ^(x). (Strictly speaking, the real function [H^ (a))H^(a))J a ^(oj) + 

(co)J^ Q ,(co>] should be considered as a whole). Thus, the effects of statistical cou- 
pling on the measured displacement spectrum will depend on the damping and natural fre- 
quencies of the modes, the excitation cross-power spectral density, and the location at which 
the measurements are made. 


As an illustration of the effects of the cross terms, the response of a simply supported 
plate to turbulent boundary layer excitation has been calculated using the methods of 
reference 6. The plate is assumed to be aluminum with dimensions of 30.48 by 1 7.78 by 
0. 1 02 cm ( 1 2 by 7 by 0.04 in.) and a damping ratio that is the same for all modes. The panel 
is exposed to typical Mach 2.0 flow. Displacement spectra were calculated both with and with- 
out the cross terms. The contributions due to statistical coupling at the eight lowest natural 
frequencies of the panel are shown in figure 2. For a damping factor of 0.04, the contribu- 
tions due to statistical coupling vary from about 10 - ^ to 10”^ times the total power spectral 
density (including cross terms). In all cases, the contribution due to the cross terms decreases 
with damping ratio. 

The influence of the cross terms on the measured damping can be estimated by compari- 
son of the calculated spectra. It is found that the error is approximately proportional to the 
cross-term contribution shown on the ordinate scale of figure 2. For the example under con- 
sideration and a damping factor less than 0.02, the error in the estimated damping factor will 
be less than 7%. 

It is interesting to note that the contributions due to statistical coupling, shown in 
figure 2, are highest for modes (1,2), (2,2), (4.1), and (2,3). Wilby (ref. 6) has shown that, for 
the boundary layer model used in the above calculations, 

J a p = 0, for n + s odd 
#0, for n + s even 

where mode orders a = (m,n) and P = (r,s). Adjacent modes that satisfy the conditions that 
(n + s) be even will be expected to show the effects of statistical coupling most clearly. For 
the modes identified above, mode pairs (1,2), (2,2) and (4,1), (2,3) satisfy the condition (n + s) 


even. 



4. 1 .4 Background Noise and Off-Resonant Contributions 

The measured displacement power spectral density function can be influenced by 
several undesirable effects, and the signals causing these effects can be referred to as noise. 

The noise may be in the form of additional excitation forces or electronic noise in the measur- 
ing equipment. Off-resonant contributions from neighboring modes also come within this 
definition of noise and can be treated in a similar manner. 

Consider first the effects caused by a second excitation that is additional to the turbu- 
lent boundary layer pressure field. For the sake of simplicity, assume that the system is a 
single-degree-of-freedom dashpot and the input and output signals are related as shown in 
figure lb. The two excitations are denoted by pj(t) and p 2 (t), where p 2 (t) represents the 
noise signal. 

Then 


S w (co) = |H a M| 2 P(w)P*(co), 


where 


P(w) = JT[ Pl (t) + p 2 (t)]e- i0Jt dt 
= PjCcd) + P 2 (cu). 


Thus, 


S w (to) = |H a M| 2 [P,(cu) + P 2 (w>] lP,(ou) + P 2 (cu)l * 
If pj(t) and P 2 (t) are uncorrelated, 

P j (o>)P 2 *(cu) = 0 = P , *(tu)P 2 (cu) 


and 


S w (cu) = |H a (co)| 2 1 


S p2 (co) 

S p ,(«) 


Sp i (cu) 


= K'(co)|H a (w)| 2 S pl (cu), 


( 21 ) 


where Spj(cu) and Sp 2 (co) are the power spectral density functions associated with p j(t) and 
P 2 (t), respectively. If Spj(cc) and Sp 2 (co) are slowly varying functions of co, then even in the 
presence of noise in the excitation, the displacement power spectral density function is 
approximately proportional to the squares of the frequency response function. Under these 
conditions, noise in the excitation introduces only negligible errors in the damping measured 
using the response power spectral density method. 



Extending the above analysis to the case of a single mode of a rectangular panel, 
equation (21) becomes 

S w (x,<4 a = ^ a 2 (x)|H a (cu)| 2 ^^ a (x')^(x")[S pl ^x",o;) + S p2 (x',x w ,w)] dx'dx", (22) 


where S p ](x',x",to) and S p 2 (x',x'',to) represent the cross-power spectral density functions of 
the two signals, which are again assumed to be mutually uncorrelated. Thus, if damping 
estimates are now made directly from S w (x,co) a , the displacement power spectral density for 
mode a, it is required that both S p j(x',x",cj) and S p 2 (x',x",o>) be slowly varying with fre- 
quency. (In fact, the requirement applies only to the real part of each cross-power spectral 
density function.) 


When noise is present in the response signal also, the measured response, for a single- 
degree-of-freedom system, has the form 


S w (co)q; 



S p2 (w) 

S pl (co) 


„(w) 

W r Pl 


= j K' (co) | H a (cu) | 2 + K” (u>) j S pl M, 


(23) 


where the noise signal n(t) is uncorrelated with the response w(t). Thus, noise in the response 
introduces an additional term K"(cu) in the total response power spectral density function. 
The displacement power spectrum can be used for an accurate estimate of the damping ratio 
only when K”(w) can be neglected or when corrections can be made to the data. 


The effects of contributions from the off-resonant modes will be similar to those of 
background noise in the response. If the effects of statistical coupling are neglected, the dis- 
placement power spectral density function for two adjacent modes, with no other noise in 
the system, is obtained from equation (17). 


S w (x,cc) = ^A 2 S p (cu)j |H a (co)| 2 ^ 2 (x)J aa (oj) + |H^(co)| 2 ^ 2 (x)J^ (w)| 
= j K a |H a (co)| 2 + KplHpMI 2 ) S p M. 


(24) 


Comparing equations (23) and (24) the term for mode P in the two-mode system has 
the same effect as the noise term K w (tu) in equation (23). Therefore, these two sources of 
error can be investigated using the following representation for the displacement power 
spectral density function: 


S w (x,t4 


1 +a(l%) 7? ] |H a (cu)| 2 K (w)S p (w), 


(25) 


where a denotes the ratio of background noise power spectral density to the true response 
power spectral density ,at resonance frequency co a . For V= 0, the background noise has a 



constant power spectral density. Other frequency dependencies for the noise can be obtained 
by assigning appropriate values of rj. For example, in the case of two modes, the contribution 
of the off-resonant mode 0 can be approximated by assigning positive and negative values to 
V for fp > f a and f^ < f a , respectively. 

The errors that the response power spectral density method can introduce into measure- 
ments of the damping ratio when background noise is present have been estimated for values 
of rj equal to 0, +2, and ±4 and for different levels of background noise. These errors are 

shown in figure 3, as a percentage of the true damping factor, with? a ,rj, and a 
as parameters. For damping factor < 0.04, the error is almost independent of and rj 
and can be treated as a function of a only. 

Figure 3 cannot be used directly to correct measured damping factors in practice. A 
typical response spectrum will contain a series of peaks and valleys. In such a case, the ratio 
of power spectral density in the valley to the resonant power spectral density of the adjacent 
peak can be used to obtain the level of background noise. A typical chart that can be used to 
correct the measured damping factor, when is less than 0.04, is contained in figure 4. 

The above procedure will not give good estimates of the error when the two natural 
frequencies f a and f^ are very close together. Under such conditions the contribution from 
the off-resonant mode will not be represented very accurately by a function of the form co^. 

It is interesting to note that, in the present case, where the vibration and noise are both 
random, mutually uncorrelated signals, the power spectral components can be added together 
without phase problems. This is an advantage over the deterministic case where the phase 
difference between resonant and off-resonant contributions is important and changes rapidly 
in the neighborhood of the resonant frequency f a . 


4.1.5 Filter Bandwidth 

One of the critical problems in the use of the response power spectral density method 
is the finite bandwidth of the filter system, whether it be analog or digital. The finite band- 
width of the filter will cause loss of resolution in the spectra with associated broadening of 
the measured bandwidth of the resonant peak and overestimation of the damping-factor. 
Filter bandwidth effects have been thoroughly investigated by Forlifer (ref. 1 1) for viscous 
damping and by Wilby (ref. 6) for hysteretic damping. Therefore, detailed analysis of the 
subject is not required here and only relevant results are presented. 

Forlifer assumes rectangular filter characteristics so that the filter transfer function 
<J>(f) has the form 

z\fp A fp 

*(0=l fp-- 4 ^ 1 <f<f F + 3 1 

= 0 


elsewhere. 



The center frequency of the filter is fp and the bandwidth is Afp. Now, if Afjyj is the 
measured half-power bandwidth of the resonant peak, which has an actual bandwidth Af^, 
then Forlifer (ref. 1 1) shows that Afp, Afj^, and Af^ are related by the equation 


tan 


-1 




This reduces to the simpler form 



(26) 


(27) 


if the trivial solution Afp/Afy^ = 0 is excluded. 


For practical application, equation (27) can be rewritten in a more convenient form 
that relates the filter bandwidth to the measured bandwidth. Then 


( ? m) = ( a, m) = (' ~(Af M /Af F ) 2 )’ <28) 

where is the ratio of the true damping to the measured value. The curve for equation 

(28) is plotted in figure 5 and can be used to correct measured damping factors. As an example 
from the curve, the error in the uncorrected measured damping factor will be less than 14% if 
the filter bandwidth is less than half of the measured resonance half-power bandwidth. 

From figure 5, it is obvious that the narrower the filter bandwidth, the smaller the error 
in the damping factor. However, there are other constraints on the filter bandwidth, the 
most important of which is the statistical confidence in the data. The statistical reliability of 
a spectral measurement can be defined in terms of the normalized standard error e = 1 fj AfpT 
where T is the total sample length. Thus, for a constant value of e , or statistical reliability, the 
total sample length must increase when the filter bandwidth is decreased. In practice, the true 
bandwidth of the resonant peak will not be known until measurements have been made. Thus, 
some form of iteration will be necessary in the selection of the filter bandwidth. The filter 
bandwidth will also have to be a compromise between the resolution and statistical confidence 
requirements. 

Additional discussion of the sample length will be presented in section 7.0, along with 
other requirements on data analysis procedures. However, it can be noted at this stage that 
the finite filter bandwidth problem is analogous to the truncation errors associated with 
some of the other methods. 



4.1.6 Drift in Natural Frequency 


A vibrating thin panel will generate heat due to internal dissipation of energy and will 
cause temperature changes in the panel. The panel is generally surrounded by relatively thick 
supports that will experience much smaller temperature changes. The differential tempera- 
ture will cause thermal stress in the panel, and the resonance frequencies of the panel will 
change accordingly. Such a drift in natural frequency was noted by Wilby (ref. 1 2) and was 
also observed in the demonstration experiments of the present investigation. An exhaustive 
investigation of such a phenomenon should include considerations of the thermal equilibrium 
and thermal stresses of the panel and their effects on the vibration characteristics of the panel. 
Such a complete investigation is beyond the scope of the present project and, instead, an 
intuitive, workable model is presented here for investigation of the effect of drift of the 
resonance frequency on the measurement of damping. 

A linear system with constant natural frequency will have a frequency response function 
H a (u>) that is independent of time. Drift in the natural frequency cu a will give a time-varying 
frequency response function H^u^t). The effect of such a drift can be easily observed 
experimentally using discrete frequency excitation. When the panel is excited by a constant- 
frequency pure tone, the drift in natural frequency will give a change in vibration amplitude 
at the excitation frequency. When random excitation is used, the response amplitude will 
have random phase that will mask the change in amplitude due to the drift in natural fre- 
quency. Thus, the presence of drift in natural frequency will not be readily detected when 
using random excitation. For a time-invariant system (as shown earlier), the response power 
spectral density function will be proportional to |H a (co)p, whereas, for a system with time- 
varying frequency response function H a (co,t), the response power spectral density function 
will be propor tional to the time-averaged (or mean) square of the frequency response func- 
tion |H a (cj)|2. 


|H 0 (co)| 2 =^/ |H a (oo,t)| 2 dt, 

1 o 

where T is the data sample length. 

Let denote the equilibrium (time-independent) resonance frequency of the system. 
Then the drift in natural frequency can be represented as 

w a (t) = w' a +w'(t), 


where 


a/(t) << for all t. 


Substituting co a (t) in the frequency response function and neglecting terms of the order 
(cu'(t)) 2 and higher, 


|H a (cu,t)h 


M a 2 [A + Bco'(t)] 


(29) 



where 


A = (w ' a 2 - w 2 ) 2 + 4 f a 2 oj^ 2 u> 2 

B = 4co^[a; 2 (2r a 2 - 1)+ w' a 2 ]. 

As an example, |H a (co,t)| 2 has been estimated for the following form of co'(t) 

cj'(t) = cj 0 'e - ^^ for t ^ 0 

= 0 for . t< 0. 


Assuming small damping factors, 


|H a (w,t)| 2 = |H a (co)| 2 (l +ej), 

where |H a (co)| 2 = l/(M a 2 A) = square of the amplitude of the equilibrium frequency 
response function and error =(l/kT)log e [(A + B co 0 e"kT)/(.A + Bco 0 )] . 

Substituting for B and A in the expression for ej, it can be seen that the error depends 
on kT, nondimensional frequency drift co'/gj^ , and the damping factor The error ej is such 
that, in general, drift in the natural frequency will produce an average frequency response 
function with a bandwidth that is wider than that for the corresponding time-invariant fre- 
quency response function. 

In practice, the error can be reduced by minimizing frequency drift using two approaches. 
First, by keeping panel deflection at the smallest possible level, the heat generated in the panel 
and the frequency drift can be kept to a minimum. Second, the thermal time constant of the 
panel can be reduced by blowing low-velocity air over the panel. This is equivalent to increas- 
ing kT in the above discussion. 


4.1.7 N-Point Measurements 

In general, the power spectral density method of determining damping factors uses a 
single measuring location. The problem of off-resonant vibration in adjacent modes, and 
background noise, could be alleviated by the use of several measuring transducers that record 
the response simultaneously. 

The displacement cross-power spectral density of an L-degree-of-freedom system can 
be written as 


L L 

s y = ^ A 0* .\b . p 
°w a — i mi mu 1 mn 

n=l m=l 


(i j = 1,2, . , . L), 


(30) 
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where 


and I mn is defined as in equation (18). 

In equation (30) S w y= S w (x*,x J ,w) and ^ n j = ^(x-*). 

Now, Xjj is a unitary matrix such that 

L L .. # 

7 -> ^ X pi S w Ux jq = r? pp 5 pq, 

j=l i=l 

where 

5 pq =1 for p = q 

= 0 forp^q. 

Multiplying the right side of equation (30) by the appropriate unitary matrix, 

L L 

^ppSpq “E S ^pn ^nm/Jmq, 

n=l m=l 


( 31 ) 


(32) 


(33) 


where 


Ppn 


i=l 


pi 


(34) 


For a single-degree-of-freedom system, 

n _ Y c £ 

1 nm A a°na ma 


(35) 


and equation (33) reduces to 

^pp^pq ~ ^a^pa^aq or *?pp “ ^a^aa^aa^ 06 ) 

For a general case where r mn can be represented by sum of L modes (corresponding 
to L degrees of freedom). 


L 

^mn ~ 2 ^i ^ni ^mi 
i=l 


( 37 ) 



and 


^pp^pq 2 ^pi^iq 
i=l 


If the response is now measured at only N different positions, the cross-power spectral 
density functions can be represented by an [N x N] matrix, i.e., the measuring system will 
allow only N degrees of freedom and there are only N distinct values of hpp’s. Thus, from 
equations (35), (36), and (37), the following conclusions can be drawn: 

a) For L = 1 , all the power will be concentrated in only one mode. 

b) For N, one can separate the power in each mode. 

c) For L > N, the scheme will fail to separate contributions from different modes. 

It should be noted that power r 7 pp represents the sum of the response power spectral 
densities in a particular mode at the N measuring locations. In practice, in a given frequency 
range, there are only a few modes with significant contributions. Thus, the condition L < N 
can be effectively satisfied. In practice, with the advances in digital data processing, the 
response power spectral density matrix can be calculated readily. The matrix diagonalization 
can be achieved easily with the use of a computer. 

The method appears to have certain advantages for resolving close natural frequencies 
and damping factors in the presence of high background noise. For all other applications, the 
method is considered to be unduly complicated when compared with the single-point 
measurement of the power spectral density method. 


4.1.8 Summary 

The basic characteristics of the displacement power spectral density method for damp- 
ing measurement have been discussed with particular reference to a panel exposed to a mov- 
ing airstream. In principle, the method is straightforward, requires no excitation other than 
the turbulent airflow, and has no instrumentation requirements other than those normally 
associated with the measurement of panel response to turbulent airflow. However, there are 
a series of problems in practice— statistical coupling, off-resonant vibration, background noise, 
filter bandwidth, and drift in the natural frequency— that have to be evaluated with respect 
to other methods. 

The influences of these problem areas on the accuracy of the damping measurements 
have been discussed, and methods of alleviating the problems have been proposed. Suitable 
positioning of the measuring transducer can emphasize the vibration in the mode of interest, 
and errors caused by background noise and filter bandwidth can be corrected by the use of 
appropriate correction charts. For optimum location of the transducer, some initial knowledge 
of the panel mode shapes is required. This is not difficult for laboratory structures but may 
pose problems for large, complicated structures of typical aerospace construction. 
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In general use, only the response power spectral density is measured. This is adequate if 
the excitation power spectral density is constant (or approximately constant) in a frequency 
range close to the natural frequency of interest, but measurements should be made of the 
pressure field to confirm the assumption. 

The use of a multipoint measuring array has been proposed as a way of reducing the 
effects of background noise and close natural frequencies. However, the use of such arrays 
involves considerable analysis problems that detract from the value of the method for general 
application. 

One problem, the drift of panel natural frequencies, is really a problem of the test speci- 
men rather than the analysis method. The drift would introduce measurement errors in all the 
test methods and the only satisfactory solution lies in the correct design of the experiment. 


4.2 RESPONSE AUTOCORRELATION FUNCTION 


4.2.1 Introduction 

The response autocorrelation method uses the decay rate of the autocorrelation function 
as a measure of the damping of the system. Thus, the approach is similar to the deterministic 
method where the damping is estimated from the reverberant decay of the vibration when the 
excitation is removed. For the ideal case of a single-degree-of-freedom system, the method 
has obvious practical value, but for multimodal systems, such as a panel, the method has 
several limitations. These limitations will be discussed in the analysis where it will be assumed 
that the vibration of the panel results from exposure to the airflow and that no secondary 
excitation is intentionally introduced. 


4.2.2 Basic Analysis 

The basic theoretical analysis for a rectangular panel exposed to random excitation, 
particularly that due to a turbulent boundary layer, has been discussed in detail by Rechtien 
(ref. 13) and Wilby (ref. 6). Consequently, the present analysis will provide only an outline 
of the method. 

The displacement autocorrelation function for the panel has been shown in equation (14) 

to be 




The integration with respect to to can be performed by contour integration, provided 
that Sp(x',x",oo) is of order less than co 2 . When t >0, 


MS p (x ',x" ,co)e iwT dco = 



where 


+ S a * cos(co d r 


+ d a 



(39) 


S(x - S a r (x ,x”, W (j + 

s a = S a *(x ',x ", a> d + iw a ? a ) 


S p (x',x ", w d + iw a f a ) = S a r - iS a ‘ 


D a/J D a@ + iD aj3 


= (w a 2 - oy 2 ) - 2w a r a (w a ? a + ?p) + 2iw d (w a r a + 


Sin ft 


cos0 =.. D si. 

IF\ Vno I r\ I 


w m 


a/3 1 


|D ap | 


and co d = Ji - is ^g na t U ral frequency for damped free vibrations. 


To advance further, assume that the damping is small so that oo a $ a « cod anc * 
S p (x',x",co d + i^a^a) can be replaced by S p (x',x",co d ). Then S a r becomes C p (x',x",to d ) 
and S a ' becomes Qp(x',x",co d ), where 

S p (x',x",oo) = C p (x',x",co) -iQ p (x',x",co). 

Also, since, by definition, the autocorrelation function is symmetrical about r = 0, the solu- 
tion to equation (38) can be extended to include all r if t is replaced by | r| in equation (39). 
Thus, 


R w (£,r) * a (x) ! c p(5'.r.«d)^“dW + %) 


+ Qp(x',x",co d )cos(co d |r|+ 0 a(3 )jdx'dx" 


for all r. 



It has been shown (ref. 6) that 


^ a (x")Qp(x',x" ) cj)dx'dx" = 0, 


so that 


R W ^’ T > = ^ 


7re 




a 2M a 2 co d oJ a 2 ? a 


,co d > sin (co d |x| + 0 aa )dx'dx" 


*<**¥*> i 

!Cp(x',x",a> d )sin(cJdkl + 6 a n) + Q p (x',x",(o d )cos(a) d |r|+ 0 a/3 )|dx'dx" 


(40) 


for all r. 


Equation (40) shows that the response autocorrelation function for the panel displace- 
ment is a sum of contributions from all normal modes plus additional cross terms associated 
with the statistical coupling of the modes. Some simplification has to be achieved before the 
method has a reasonably practical form. 

From (ref. 6), when the natural frequencies are well separated and the damping is small 
(f << 1), the cross terms can be neglected. When the natural frequencies are close together, 
the conditions for small cross terms are more difficult to establish but are similar to those 
associated with the displacement power spectral density function. 


Neglecting cross terms equation (40) becomes 


R w (x,r) 


7re CO « r <* |r| 

“ 2M a 2 co a 3 r a d - f a 2 > 1/2 


^ a 2 (x)sm(w d |r| + V) 




(41) 


From equation (41), the autocorrelation function is a sum of the contributions from the 
different modes and it is difficult to estimate the damping coefficient unless the vibration due 
to one mode can be isolated. This can be achieved to a certain extent if the natural frequencies 
are widely separated and the vibration in the vicinity of one natural frequency can be effec- 
tively isolated by filtering. The integral term in equation (41) is the joint acceptance term in 
equation (17) for the displacement power spectral density function. 



If the vibration for one mode can be isolated, then the autocorrelation function becomes 
an exponentially decaying sine function where the decay rate is proportional to the damping 
and the sine function has a nonzero phase angle. For small damping d aa »n/2. 

From equation (41), the autocorrelation function for a single-degree-of-freedom system, 
such as a dashpot, is 


R w (t) = 


+ 0 aa ) 

2M a 2 co a 3 f a a - r a 2 ) 1/2 


S p (co d ). 


(42) 


4.2.3 Multimodal System With Close Natural Frequencies 

The autocorrelation function R w (x,r) for a multimodal system has been obtained in 
approximate form in equation (41) where the cross terms (a^0) have been neglected. When 
a single mode can be considered alone, the autocorrelation function can be used to calculate 
the modal damping, but when there are two or more modes with close natural frequencies, 
modal isolation is not possible. 

Consider the case of two modes of a panel where the natural frequencies are closely 
spaced. Using the approximate form for the displacement autocorrelation function (equation 
(41)), the autocorrelation for the two modes can be written as: 

R w (x> r ) = A a sin ^da M + d a0i ) + A0sin(w d 0|r| + dpp), (43) 


where 


Aq, 


i w -0J a?al r l 
(x)ne 


2M a z cVr a (l - r a z ) 




2 Vi Ja j A 


and 



Similar definitions hold for A^ and 

To illustrate the relative importance of coefficients A a and Ag, assume that 6 aa = dpp = 6 
for closely spaced modes since it is unlikely that the damping of the two modes will differ by 
a large amount (d aa = tan~l|(l - fa 2 )^/^)- Then, for t> 0, equation (43) becomes 

R w (x,t) ss A a sin(co da T+ 0) + A^sin(co dj3 t+ 0) 


where A^and A^ are always positive. 



Letting c opp = w da + co', then 

o') 

R w (x,r) = (A a z + Ap Z + 2A a A^cosc o'r) sin(co da r+ 0 + 0(r)), (44) 

where tan0(r) = (A^sinw'r)/(A a + A^cosgj't). 

Thus, 0(r) acts as a varying phase angle and 2A a AjgCOSGJ 'r causes beating of the ampli- 
tude of the autocorrelation function, varying from A a + Aq to A a - Ap. The maximum error 
in the autocorrelation amplitude will be ±A^. These results are similar to those of White (ref. 
14) for the impulse response function. 

To demonstrate the effect of close natural frequencies on the displacement autocorrela- 
tion function, equation (43) has been evaluated for two modal pairs of a 30.48- by 17.78- 
by 0.102-cm (12- by 7- by 0.04-in.) simply supported panel exposed to Mach 2.0 airflow and 
a boundary layer displacement thickness of 0.589 cm (0.232 in.). The method of reference 6 
was used to calculate the terms A a and Ap. In one example (fig. 6), the mode pair is a= (3,1), 
(3 =0,2) and the damping is assumed to be = 0.007 = For the second example (fig. 7), 
a = (2,2), (3 = (4,1), and the damping is again = 0.007 = $p. The data in figures 6 and 7 are 
plotted in terms of the normalized autocorrelation function or autocorrelation coefficient, 
and the autocorrelation curves are compared with the corresponding curves for mode a alone. 

Figure 6 clearly shows the beating between the contributions from the two modes, and 
experimental data of this type would indicate the presence of more than one mode. However, 
in figure 7, for the time-delay range shown, the beating phenomenon is not evident, and 
experimental data in this form could be interpreted as being associated with a single heavily 
damped mode. Thus, it is essential that, for proper evaluation of the data, the maximum time 
delay t in the autocorrelation function must be sufficiently large to indicate the presence of 
more than one mode. To identify one period in the beat signal, T m must be greater than 
2jr/to'. In figure 6 , 27t/cj' = 12.45 milliseconds, r m is 17.3 milliseconds, and the beat can be 
identified. In figure 7, 27 t/cj' = 22.57 milliseconds and T m is 12.6 milliseconds, so the beat 
cannot be identified. 

In practice, the dominant factors in the coefficients A a and A^ will probably be the 
mode shape functions \p a (\) and 0o(x), and choice of measuring position x will have an 
important influence on the form of the autocorrelation function. For the examples in 
figures 6 and 7, a location xj = 5.08 cm (2 in.), X3 = 4.45 cm (1.75 in.) was chosen. 


4.2.4 Band-Limited Signal 

For the multimodal system, it has been shown above that the presence of more than one 
natural frequency has a distorting effect on the single-mode autocorrelation function. Even 
in the case of two modes, the distortion can cause large errors in the damping estimates. 

In practice, the effects of neighboring natural frequencies can be minimized by the use 
of filters, although this method has very limited application. The filters exclude vibration at 
frequencies away from the natural frequency of interest and should have steep cutoff 



characteristics. However, since the filters cannot exclude off-resonant vibration contributions 
at the natural frequency under investigation, they can be used successfully only if the ratio of 
resonant to off-resonant vibration is large. The use of filters may introduce adverse effects on 
the autocorrelation function if the bandwidth is too narrow. Thus, an optimum filter band- 
width has to be selected. 


The effect of filter bandwidth on the autocorrelation function of the mode of interest 
can be investigated using a single-degree-of-freedom system as the model. The narrowband 
vibration signal could be the result of introducing a band-pass filter in the excitation or in the 
response. If the narrowband autocorrelation function is denoted by R w (r; Ac o), then 


R w (t; Aco) = I H a (oj) | 2 S p (oj)e laJT dw . 

The function 4>(a>) represents the filter, and, for a rectangular filter, 


=1 Wj < M <u> 2 


= 0 elsewhere. 


(45) 


(46) 


If it is assumed that, within the filter bandwidth, the excitation pressure power spectral 
density function has a constant value Sp, then equation (45) reduces to 

R w (r; Agj) = S /°°4»(co)|H a (co)| 2 e iwr da;. (47) 

* -OO 

Evaluation of the integral in equation (47) is discussed in appendix A. 

The problem is similar to that discussed by Kandianis (ref. 15) for the excitation- 
response cross-power spectral density method. Using the approximations of Kandianis, the 
presence of the filter modifies the autocorrelation function by the addition of terms that 
include functions of the form (cos oo 2 t)/t and (sinod 2 T)/r, which decay rapidly as r increases. 


4.2.5 Background Noise 

Background noise can occur either as an additional excitation or as electronic noise in 
the equipment used to measure the vibration. The effects of the noise signals on the response 
autocorrelation function for a single-degree-of-freedom system have been discussed by 
Kandianis (ref. 16) and the results can be incorporated in this discussion. 

Considering first the excitation, assume that a second signal p^(t) is superimposed on 
the initial signal pj(t) and that the two signals are mutually uncorrelated. Then, from equa- 
tion (16), the response autocorrelation signal is 

R w (r) = jCT h(7j)h(T 2 ) [R pl (T + Tj - r 2 ) + R p2 (T + Tj -r 2 )]dT 2 dTj. (48) 

Thus, the displacement autocorrelation function for the noise-free excitation, equation (16), 
is modified by the addition of the noise autocorrelation function to the pressure autocorrelation 



function. The change, as expected, is a parallel of the effect demonstrated in equation (21) 
for the displacement power spectral density. 

If the signals pj(t) and P 2 (t) represent white noise, 

R pl (r) =27rS pl 5(r), 

R p2 (r) = 27rS p 25(T), 


and 


R w (t) = 2tt(S p j + S p2 )Pjp h(rj )h(T 2 )5(r + r, - r 2 )dr 1 dr 2 . (49) 

In equation (49) S p j is replaced by [S p j + S p 2l , which has the effect of simply magnifying 
the autocorrelation function. 

As an alternative to equation (48), the autocorrelation function can be written in the form 

R w (r) = J[jH a ( w )l 2 [S pl M + S p 2 (w)] e iwT dw, (50) 

where there are now no restrictions on the excitation power spectral density functions, and 
[S p j(cu) + S p 2(co)] replaces the term S p j(w) found in the noise-free equations. 

Referring now to the multimodal system associated with a panel and assuming that the 
statistical coupling cross terms can be neglected as in equation (41), the effect of the second 
(noise) excitation is determined by the real part of the cross-power spectral density function. 
Using equation (41) as the reference, the integral 

/ / tf' a (x')\l' a (i") C p (x', x",a>)dx'dx", 

l\ l\ 

associated with the noise-free case, becomes 

J[ f A <y>0* a (x")[C pl (x', x-.w) + C p2 (x', x",co)]dx'dx« 

when noise is present in the excitation. Here the effect of the noise signal depends not only 
on the power spectral density, but also on the spatial correlation in relation to the mode 
shapes of the panel. The function C p 2 (x', x",w) represents the real part of the complex cross- 
power spectral density function associated with the noise signal. 

When noise is present in the response system, there will be no influence on the panel 
vibration, and the noise will appear as an additional term in the autocorrelation function. If 
R n (r) is the autocorrelation function associated with the noise in the response system, and if 
the autocorrelation function for the vibration is denoted by R W '( T )» then the measured auto- 
correlation function will be 


R w (t) = R w '(t) + R n (r). 


(51) 
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If the noise signal n(t) is broadband white noise with power spectral density S n , then 


R n (r) = 27rS n 5(r), (52) 

where 8(r) is the delta function. Therefore, the measured autocorrelation function differs 
from the noise-free correlation only when r = 0, at which time delay the noise introduces a 
delta function spike. 


In practice, it is more likely that the noise will resemble band-pass noise rather than 
wide-band noise. The autocorrelation function then becomes 


sin c ojT sin ojiT 

R nW = 2S n |— ,ri- 


(53) 


where an d are the upper and lower cutoff frequencies, respectively, for the band-pass 
filter. If noise of the type represented by equation (53) is present in the measurements, it 
will be difficult to make suitable corrections to the data. Thus, if R n (r) is a significant fraction 
of R w (r), the damping measurements will show appreciable errors. 


4.2.6 Summary 

The response autocorrelation method involves analysis of the vibration in the time do- 
main, and this results in additional complications when compared with frequency domain 
analyses. Analyses of the basic characteristics of the method show that severe problems will 
occur in practice unless the natural frequencies of interest are well separated. In such cases, 
the vibration in the neighborhood of a given natural frequency is due essentially to only one 
mode. Unfortunately, these very restrictive conditions are satisfied by very few modes of a 
structure, even a simple panel. 

In an attempt to isolate the vibration of one mode, the use of band-pass filters has been 
suggested. Here a compromise has to be reached so that the filters have a bandwidth wide 
enough not to distort the results for the mode of interest, yet narrow enough to exclude 
vibration from other modes. Irrespective of filter characteristics, they will not exclude off- 
resonant vibrations within the band-pass frequency range. Thus, the filters are of little use 
unless the resonant vibration is large with respect to the off-resonant vibration. 

Background noise, either in the excitation or the response, can cause errors in the damp- 
ing measurements. This is particularly true if the noise is in the response signal. However, if 
the noise in the excitation has a constant power spectral density within the frequency range 
of interest, the associated errors in the measured damping will be negligible. 



4.3 EXCITATION-RESPONSE CROSS-POWER SPECTRAL DENSITY 


4.3.1 Introduction 

The excitation-response cross-power spectral density method was proposed by Clarkson 
and Mercer (ref. 17) and is the stochastic equivalent of the response amplitude phase method 
(see Kennedy-Pancu, ref. 4) for discrete frequency excitation. Phase information, which is 
not used in the response power spectral density or autocorrelation methods, is used in the 
cross-power spectral density approach. Therefore, at least in principle, the method is 
capable of separating contributions from off-resonant modes and background noise. However, 
in practice, several limitations are imposed on the method, and these limitations are discussed 
in this section. 

When studying structural response in the presence of an airflow, a severe limitation 
arises because of the inability to measure the excitation and response fields simultaneously. 
Introduction of a microphone would disturb the flow field or require modifications to the 
structure. Thus, it is necessary to provide a second excitation that can be used as a reference. 
The pressure field due to the airflow then takes on the role of unwanted noise as far as the 
damping measurement technique is concerned, and the problem of signal-to-noise ratio 
becomes very important. 


4.3.2 Basic Analysis 

The basic analysis for the method will be performed under the assumption that there is 
only one excitation present, and this excitation will be used as the reference signal. Thus, 
the system is identical to that shown in figure la where the excitation signal is p(x',t) and 
the response of the system is denoted by w(x ,t). It is assumed in the analysis that the excita- 
tion is applied at a single point rather than being distributed over the panel, but this assump- 
tion is not necessary in general. The cross-power spectral density function relating the excita- 
tion and response can be represented by the equation: 

S w , p (x, x', to) = W(x, to) P*(x ',«), (54) 

where W(x, to) and P(x',co) are the complex Fourier amplitudes of w(x, t) and p(x', t) 
respectively, and the asterisk denotes the complex conjugate. From equations (8), (9), (10), 
and (11), 


W(x,co) = 52 *// (x) H (to) f i// (x) P(x,to) dx. 
oc « A u 

Substituting for W(x,to) in equation (54), 

S w , p (x, x',w) = £H a (co) * a (x) J^S p (co) tf a (x)5(x' - x) dx 
= EH a Mi// a (x)^(x')Sp(co), 

where Sp(co) is the power spectral density of the pressure signal at x'. 



The right side of equation (55) indicates that the cross-power spectral density function, 
which is complex, depends on the frequency response function H a (co), the mode shapes 
\p a (x) of the panel, and the excitation force. The two mode-shape terms will give a multi- 
plying factor that determines the efficiency with which the vibration is being excited and 
measured. For example, if the point of excitation and the location of the measuring probe 
are at antinodes of the mode being measured, then ^ a (x') = ^ a (x) = 1 and the measured 
cross-power spectral density function will be a maximum for a given excitation force. 

Conversely, if the exciter or measuring transducer is located at a node of a given mode, 
then no vibration will be recorded in that mode. These properties are valuable in practice for 
augmenting or suppressing the response in a particular mode. 

If an excitation force with constant spectrum level is used, then Sp(co) will be independent 
of co and the excitation-response cross-power spectral density function will be directly pro- 
portional to the sum of the frequency response functions of all modes. Comparing equation 
(55) with equation (20), note that the cross-power spectral density function does not contain 
the double summation present in equation (20) for the response power spectral density func- 
tion. Cross terms that arise from the statistical coupling in the power spectral density function 
are absent in the excitation-response cross-power spectral density. However, the frequency- 
dependent term H^co) will still impose filter bandwidth requirements if adequate resolution 
is to be achieved in the data reduction of the cross-power spectra. 

At some stages in the subsequent analysis, it will be convenient to refer to the non- 
dimensional cross-power spectral density function T w p(x, x',co), which is defined as 


n , , . s (x,x,w) 

r w n(x.x > W ) = ■ lP I7T 

W ’ P [S w (x,co)-S ( CO )]*/ 2 


(56) 


The function T^, n (x,x',co) is a measure of the coherence between the excitation and response. 

w,p — — 

For a single-mode system, the dimensional and nondimensional cross-power spectral 
density functions for a single excitation and the associated response are 


S wp (x,x',co) Q; = H a (co) \J/ a (x) (// a (x')S p (to) 


(57) 


and 


r W;P (x,x',w) a = H a (w)^ a (x),// a (x') 


S p M 
S\y(x> 


T 1/2 


( 58 ) 


Since S w (x,co) a , the displacement power spectral density in mode a, is proportional to 
|H a (cc)| 2 , then r w>p (x,x', 0 L>) is proportional to H a (oj)/|H a (eo)|. 



4.3.3 Background Noise 


As stated earlier, one important problem when applying the excitation-response cross- 
power spectral density method to a panel exposed to turbulent airflow is the difficulty in 
measuring the excitation. To overcome this problem, the introduction of a second excitation 
has been proposed (refs. 6 and 18) as the reference signal. The system then resembles that 
shown in figure lb where the boundary layer pressure field p j (x',t) is regarded as undesirable 
noise and the signal p 2 (x",t) becomes the reference signal. 

The modified cross-power spectral density function can be determined following the 
approach of section 4.3.2, except that the panel total displacement now depends on two 
excitations, only one of which (p 2 ) is being measured. If pj(x',t) and p 2 (x",t) are assumed 
to be statistically independent, the nondimensional cross-power spectral density function 
for a single mode becomes 

[S d2 ^>] 1/2 »«(«) 

r w n -)(x,x",oS) = & — x— ^ =rn^ • (59) 

W ’ p2 ■ - “ [A 2 S p] (w)J aa (w) /4+ </z a 2 (x") S p2 (co)] !/2 * a (x)|H a (w)| 

When [A 2 S p jMJ aa (co)/4] is much smaller than [i// a 2 (x")S p2 (a))] , which means physically 
that the panel vibration is due mainly to the introduced excitation p 2 (x",t), then 
T w p2(£’5 w > tJ ) i n equation (59) tends to the noise-free form in equation (58). 

If noise n(t) is introduced into the response, as indicated in figure lc, and n(t) is 
assumed to be uncorrelated with the vibration w(x,t) and the excitation p 2 (x",t), then the 
measured nondimensional cross-power spectral density function for mode a will be given by 
the expression 


lMx)* a (x") [S 2 (co)] !/ 2 H (oo) 

T w p2 (x,x",co) l a = T- x x f ^ j-Z TTT) ‘ (60) 

,P j|H a M| 2 ^ a 2 (x)[A 2 S pl (co)J aa (oo)/4 + ^ 2 (x")S p2 (co)] + S n (co) J W 2 

Equations (58), (59), and (60) all have the same numerator which, apart from the omis- 
sion of a factor [S„ 2 (co)] ^ 2 , is the dimensional cross-power spectral density function relating 
excitation p 2 (x",t) and panel displacement. The differences between the three equations occur 
in the denominator which, in the form shown, is the square root of the response power spectral 
density. Since the addition of noise in the excitation or response increases the measured 
response power spectral density, the nondimensional cross-power spectral density decreases 
as the noise increases. 

The decrease of T w p2 (x, x ",cu) a as the noise increases identifies an important problem 
area for the case of a panel excited by a moving airstream, namely, the accurate measurement 
of small correlation coefficients or functions obtained from large signals. To perform the data 
analysis with sufficient accuracy, it is necessary that the data reduction system should have 
high amplitude resolution and large dynamic range. For a given instrument system, the 
dynamic range will be constant and the maximum signal will be determined by the vibration- 
plus-noise signal. Only a small fraction of this signal will be correlated with the excitation 
reference, thereby creating a signal processing problem. 



4.3.4 Multimodal System 


The excitation-response cross-power spectral density function for a noise-free, multi- 
modal system is given in equation (55) as 

S WiP (x,!S',u) = W H a (w)S p (w). (61) 

Comparing this equation with the corresponding equation for the displacement power spectral 
density function (equation (20)), the main advantage of the cross-power spectral density can 
be identified. 

In equations (20) and (61), the value of the function at a given frequency is determined 
by a summation over all modes. For the power spectra, each modal contribution is propor- 
tional to the real function IH^fcj)!^ and is added algebraically to the other contributions. In 
contrast, the modal contributions to the cross-power spectral density are proportional to the 
complex function H a (co) and are summed vectorially. The vectorial addition is the direct 
equivalent of the Kennedy and Pancu method (ref. 4), and permits separation of the resonant 
vibration from the off-resonant/noise components. This identifies the advantage of the cross- 
power spectral density method. 

From equation (61), if Sp(cu) is essentially constant, the locus of the measured cross- 
power spectral density function will be a constant multiple of the frequency response func- 
tion obtained from the Kennedy and Pancu method with discrete frequencies. Properties of 
the frequency response functions discussed in section 3.3 will be equally applicable to the 
cross-power spectra, and resonance circles can be drawn as an aid to the measurement of modal 
damping. This technique corrects for the presence of off-resonant and noise contributions. 

The nondimensional cross-power spectral density function for the multimodal system 
can be obtained by following the procedure outlined in section 4.3.3 and using equation (61). 
The nondimensionalizing factor [Sp(co) S w (x,cu)l ^ in equation (58) will be replaced by 
[S p (cj) • g S^x,^))^] ^ , thereby reducing the value of the nondimensional cross-power 
spectral density function. 


4.3.5 Filter Bandwidth 

The fundamental quantity being measured in the cross-power spectral density method is 
the frequency response function H a (cu). Thus, one of the problems in the data reduction will 
be the loss of resolution due to the finite bandwidths of the analyzing filters. Essentially, the 
problem is similar to that encountered in the power spectral density method. In practice, 
the cross-power spectral density function can be determined in two ways that will be referred 
to as the direct and indirect methods. If digital techniques of data analysis are used, the 
direct method indicates that the Fourier transforms (or complex Fourier amplitudes) of the 
time series are computed and multiplied together to give the cross-power spectral density 
function. The indirect method first computes the excitation-response cross-correlation func- 
tion and then Fourier transforms the correlation function to obtain the cross-power spectral 
density function. Again, referring to digital techniques for data analysis, the filter bandwidth 
is relatedto the data sample length for the direct method or maximum delay time for the 
indirect method. 



Considering first the direct method, the filter bandwidth requirements can be established 
using the discussion in section 4.1.5 as a guideline. For a filter with rectangular characteristics, 
the error in the damping measurements obtained from the power spectrum is shown in figure 5 
as a function of filter bandwidth and measured half-power bandwidth. In practice, the filter 
characteristics in the direct method will be determined by the presence, or absence, of data 
smoothing but will differ from the ideal rectangular shape. Filter bandwidth is discussed further 
in section 7.0. 

The filter bandwidth problem for the indirect method has been discussed by Clarkson 
and Mercer (ref. 17) and Soovere and Clarkson (ref. 19) in terms of truncation, which is the 
equivalent time domain problem. As will be seen in section 4.4, truncation problems occur 
also in the single-sided Fourier transform method. 

In the indirect method, the theoretical excitation-displacement cross-power spectral 
density function S w p (x, x',u>) is obtained as the Fourier transform of the cross-correlation 
function R w p(x,x',fj, where 

S w ,p(x,x',cu) = J7 J[^w,p^^'’ r ) e " 1CJrdT - (62) 

However, in practice, it is not possible to measure R w p(x,x',r) for an infinitely long time 
delay, and the function must be measured for a finite range, say -r m < r< r m . The computed 
cross-power spectral density function is then 

S w,p^’S T m) = 27 /_>) R w,p(*’S , ’ T )e~ lcJ ‘ r dr, (63) 

where D(t) is a weighting function, or lag window, such that D(t) = 0 for|r| > r m . 

Alternatively, 


where 


S W)P (x,x',oj; 


T ) = f°% (W -Of) 

nv J_ oo 


S w,p^>^'’ w ') dco '> 


(64) 


4>(co -coO = 2~f°D(f)e~^ C0 _tJ ) T d r 

and is called the spectral window. 4>(co-cu0 represents the filter characteristics referred to above. 


For the case of single-point excitation, as considered in section 4.3.2, and assuming 
constant excitation spectral density, equation (64) can be written in the form 

Sw,p<x^«;im)-E K a (x,x') Haiaj'Wu)'. (65) 

The form of the lag window D(t) can be selected from several standard representations. 
Clarkson and Mercer (ref. 17) selected the rectangular lag window 


D(r) = l |r|<r m 
= 0 


elsewhere. 


( 66 ) 



As a modification, Soovere and Clarkson (ref. 19) use the Bartlett triangular window. 


D(r) = 1 -irf |T|<r m 

m 

= 0 elsewhere 



since this is more reliable, statistically, than is the rectangular window. In both cases, the 
effect of truncation is determined with respect to the parameter (l/p)(ds/df) where p is the 
radius of curvature of the cross-power spectral density curve at resonance and ds/df is the 
rate of change of arc length with frequency at resonance. Assuming small damping, « 1, 
Clarkson and Mercer (ref. 17) show that 


J_ds 4 tt 
P df 


1 - e C ° aJa ' m [ 1 + UgtS a T m + h/ 


,?/v T . 


2t 2 . 2, 


l-e' W m(|t4fc’ml 


(68) 


and a similar equation is derived by Soovere and Clarkson for the triangular lag window. Curves 
of (l/p)(ds/df) as a function of and r m can be used as a means of correcting for trunca- 
tion (or filter bandwidth) errors. Examples of such curves for the rectangular and triangular 
lag windows are contained in figure 8. The curves show that, for a given maximum time delay 
r m , the error created by the rectangular lag window is less than that due to the triangular 
window. This is expected since the equivalent bandwidth of the triangular lag window is 
larger than that of the rectangular window, for a given value of r . The advantage of the 
triangular window lies in the improved statistical reliability. 


Some simplification of the curves is possible if it assumed that the filter bandwidth Afp 
is a fraction of the resonance bandwidth Af a , i.e., 

A f p Af 0 . 

From reference 20, the effective filter bandwidth for the rectangular lag window can 
be taken as 


Then 


and 


where 


Af F-2T: 


m 


a i f t 
or a m 


= n7r/2 



*(n) = 


1 - e ~ n7r / 2 Cl + n7r/2 + n 2 7r 2 /8) 
1 _ e ‘ n7r / 2 (l + nir/2) 


(69) 
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The curve for 'l'(n) is shown in figure 9. At low values of Af a /Afp (less than 3), the 
curve for 'l'(n) indicates larger errors in the damping measurement than are predicted by 
figure 5 for the rectangular filter. 

The above discussion on truncation errors has assumed that the instrument is ideal in that 
there are no dynamic range or noise problems. In practice, there will be a dynamic range prob- 
lem when the time delay exceeds a certain value that will depend on the equipment. Any cor- 
relation information obtained for time delays in excess of the critical value will be spurious, 
and the critical time delay will determine the maximum useful time delay. 


4.3.6 Summary 

In concept, the excitation-response cross-power spectral density method has several 
highly desirable features that result from the inclusion of phase information. Using the excita- 
tion force as a reference, the response vector diagram can be drawn and the effects of back- 
ground noise and off-resonant vibration minimized. Additional benefits arise if the point of 
application of the excitation can be varied, thereby suppressing vibration in certain modes. 

When the method is applied to the case of structural response to a moving airstream, 
measurement of the excitation force becomes a significant problem. If the aerodynamic pres- 
sure field is used as a reference, a measuring transducer has to be located in the airstream or in 
the surface of the structure. In either case, the experimental environment would be changed. 
As an alternative solution, a second excitation can be introduced to act as the reference, and 
the boundary layer excitation then takes on the characteristics of unwanted noise. However, 
under these circumstances the problem of signal-to-noise ratio becomes important. The vibra- 
tion induced by the second excitation has to be significantly greater in amplitude than that 
induced by the aerodynamic pressures if reasonable accuracy is to be achieved in the damping 
measurements. To satisfy this requirement, the second excitation force must be large, with 
resultant dynamic range problems in the instrumentation. If the test panel is enclosed by a 
pressure equalizing chamber, difficulties may arise in locating the exciter within the chamber. 


4.4 RESPONSE SINGLE-SIDED FOURIER TRANSFORM METHOD 


4.4.1 Introduction 

The important characteristic of the single-sided Fourier transform method is that it uses 
phase information that is hidden in the power spectral density and autocorrelation function. 
The method was introduced by Kandianis (ref. 16) as a means of reducing the effects of noise 
interference in the analysis of transient signals. However, the theoretical background applies 
strictly to stationary signals and, for simplification, the characteristics of the signals are usually 
assumed to be those of white noise. 

In operation, the method computes the response autocorrelation function, either directly, 
or indirectly via the power spectral density function, and then performs Fourier transformation 
on half of the autocorrelation function (e.g., for r >0). The resulting single-sided Fourier 



transform is a complex function of the frequency response function H(co), and the damping 
can be estimated from the complex plane diagram in a manner similar to that developed by 
Kennedy and Pancu (ref. 4) and used in the excitation-response cross-power spectral density 
method (sec. 4.3). 

Discussion in this section will be concerned with vibration induced by random excitation 
such as the turbulent boundary layer. Application of the method to transient excitation 
techniques where a second excitation is introduced will be discussed later in section 5.0. 


4.4.2 Basic Analysis 

The basic properties of the method can be demonstrated most easily by a single-degree- 
of-freedom system. The displacement autocorrelation function for such a system has the form, 
from equation (42), 


R w< T > T£Sin( Wd |r|+S ai( )Sp( Wd ). (70) 

7MaWV' - fa) 

Then the single-sided Fourier transform is defined as 

F w ( “> * i j[ R w Me-‘" T dr. (71) 

Substituting equation (70) in equation (71), F w (ou) can be evaluated using standard 
integrals, and 
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(73) 


As expected from general principles, the real part of F w (co) is half of the displacement 
power spectral density function and is proportional to |H a (cj)|2. Thus, the real part can pro- 
vide estimates of the damping in precisely the same way as in the power spectral density 
method. However, the advantage of the method lies in the combined use of the real and 
imaginary parts in the form indicated by equation (72). 


The locus of H a (co) in the complex Argand plane has been discussed in section 3.3 where 
it was stated that H a (co) was a circle only for the case of hysteretic damping. For viscous 
damping, the frequency response function is approximately a circle in the neighborhood of 
the natural frequency when the damping is small. Similar restrictions apply to the function 
F w (co) of equation (72). 



Using the Kennedy and Pancu method (ref. 4), estimates of the damping coefficient are 
obtained from the locations of the half-power points of H^cu). This is an approximation that 
is valid for light damping and it is again valid for the function F(cu). To demonstrate the validity 
of the approximation, consider a parameter X defined by 

= ^|H a (co a )| 2 . 

Then X can be assumed to be equal to 2 (which gives the half-power law) within an accuracy 
of about r a /2 for hysteretic damping and 3? a /2 for viscous damping. In the case of the single- 
sided Fourier transform F(w), the error in X is approximately 2f a . Thus, for the single-degree- 
of-freedom system, the single-sided Fourier transform and the excitation-response cross-power 
spectral density methods are closely allied. 

The single-sided Fourier transform method can be extended to include the multimodal 
case typical of panel vibration by the use of equation (40). Substituting equation (40) in 
equation (71), 


! I W/vf/v 

H„k + -f® 



where D^, sin 0^, and cos 0 a ^ are defined in section 4.2.2 for equation (39). 

Immediately obvious in equation (74) are the cross terms (otj* (}), which are absent in the 
excitation-response cross-power spectral density function but do occur in the power spectral 
density function. The joint terms (a = /)) pose no problem since they represent a sum of terms 
of the type K a (2f a + i )H a (cj) where the parameter K a is a function of mode shape and 
measuring location and is independent of frequency co- Thus, for small damping, the locus of 
the joint terms can be analyzed in a straightforward manner using the complex Argand diagram 
and constructing resonance circles at the natural frequencies. The practical application of the 
method to a panel excited by turbulent airflow depends on the magnitude of the interference 
from the cross terms. 


When the natural frequencies are well separated, the function(cos 0 a 
^ and(sin 8 a J) |D a J is of order Thus, if « 1 , these funt 
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considered small relative to the joint term of order However, when the natural 

frequencies are close together, 6 a p and D a ^ approach the corresponding joint term values 
6 aa and D a a, and the influence of the cross terms will depend on the measuring location x 
(through i// a (x), ^(x)), and the integral terms. 


4.4.3 Background Noise 

The effects of noise on the damping estimates obtained from the single-sided Fourier 
transform method have been discussed by Kandianis (ref. 1 6) for a single-degree-of-freedom 
system. From equation (50), noise in the excitation, if it is uncorrelated with the excitation, 
is additive in terms of the power spectral density functions. Noise in the response, if it is 
uncorrelated with the response, is additive in terms of the autocorrelation function— equation 
(51). Translating these effects into the single-sided Fourier transform function of equation (72), 
then 


F w (u) = rr^tr < s Pi< w d> + Sp2<m»(2t„ + ig)+ s„m. to 

S M a 

The presence of the noise P 2 (t) in the excitation increases the coefficient of H(oo)(2f a + 
icu/co a ) but has no influence on the damping estimates. Noise n(t) in the response signal shifts 
the plot of F w (co) along the real axis. In principle, if S n (co) is constant, the accuracy of the 
damping measurements is again unaltered, but, in practice, if n(t) is large so that the dynamic 
range constraints of the instrumentation become important, accuracy may be impaired. 
Further, if S n (co) varies in the neighborhood of a natural frequency, additional errors will 
appear in the damping estimates. 


4.4.4 Truncation Errors 

Truncation errors, encountered when the excitation-response cross correlation function 
is transformed into the cross power spectral density function, are discussed in section 4.3.5 
as an equivalent of the resolution losses associated with finite bandwidth filters. The single- 
sided Fourier transform method is also susceptible to truncation errors associated with the 
transformation of the autocorrelation function. 

The truncation errors arise because it is not practical to compute the displacement auto- 
correlation function for infinitely long time delays. Thus, equation (7 1 ) for a single-degree- 
of-freedom system becomes, in practice, 

F w (cu;r m ) =-±f D(T)R w (T) e - iwT dr, (76) 


where D(t) = 0 for r > r m> 



Equation (76) is the direct equivalent of equation (63) for the cross-power spectral 
density case. Since the correlation functions in equations (63) and (76) show some differ- 
ences, equation (76) will be evaluated for the case of a rectangular lag window 


D(r) = 1 0<r<T m 

= 0 elsewhere. 



Assuming that the excitation has a constant power spectral density Sp, and that the 
system is lightly damped such that « 1 , then, from appendix B, 
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Comparing with equation (68), it is seen that equation (78) contains an extra term 
^sin^Tj^cos^tj^ in the numerator and denominator. If it can be assumed that the damp- 
ing system is so small that 


f a sinc Vm COSG Vm «*> 


then equation (78) reduces to 
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which is identical to equation (68). Thus, curves of the type shown in figures 8 and 9 can be 
used to estimate the truncation errors in the single-sided Fourier transform method, although 
it should be remembered that the condition on the damping is now, effectively, « 1 
rather than « j 


In a similar manner, the truncation errors in the single-sided Fourier transformation with 
a triangular lag window will be approximately equal to those predicted by Soovere and 
Clarkson (ref. 1 9) for the cross-power spectral density function (see fig. 8). 


4.4.5 Summary 

The single-sided Fourier transformation method operates on the response autocorrela- 
tion function in such a manner that hidden phase angle information is exposed. For the 
application considered in the preceding sections, the panel response is that due solely to the 
aerodynamic pressure fluctuations; the introduction of a second excitation will be discussed 
in section 5.0. 



The ability to use phase information places the single-sided Fourier transform in a similar 
category to the cross-power spectral density method, without the need to measure the excita- 
tion and response simultaneously. Both methods can make allowances for background noise 
and off-resonant vibration more easily than when only amplitude data are measured. In prac- 
tice, errors arising in the single-sided Fourier transform due to noise, off-resonant vibration, 
and truncation effects will be similar to those encountered in the cross-power spectral density 
method. Correction procedures developed for the latter method can be used for the single- 
sided Fourier transform although now the restrictions on small damping will be slightly more 
severe. 
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5.0 NON-STEADY-STATE METHODS 


5.1 INTRODUCTION 

The category of non-steady-state methods includes quasi-steady-state and transient 
methods. Quasi-steady-state methods are usually associated with sinusoidal excitations, which 
are slowly swept through the frequency range of interest. Narrowband random excitation 
sweeps have been proposed for structural testing (refs. 21 and 22), but their application is 
mainly that of determining response amplitude rather than the measurement of damping. 
Recently, transient methods (refs. 23 and 24) have been proposed as replacements for the 
slow sweep rates when there is a need for a significant reduction in testing time. An example 
of such a test environment is the flutter testing of an airplane, where it is difficult to achieve 
stable flight conditions in some parts of the flight envelope. Development of the transient 
method has been carried out by White (refs. 14, 25, 26, and 27) and Kandianis (refs. 16 and 
28), the latter having utilized random analysis methods to reduce the errors from noise 
interference. 

With regard to the current problem of damping measurement of panels exposed to a 
moving airflow, the use of any non-steady-state method implies the introduction of a second 
excitation at a given point of application. This has the advantage of allowing selection of a 
suitable location for the excitation and alleviates to some extent the effect of the cross terms 
of the type shown in equation (20) when a The non-steady-state methods have a second 
common characteristic because, with the exception of the narrowband random noise sweep, 
which is not applicable to the present problem, the excitation in all the above methods is 
deterministic in form. 


5.2 QUASI-STEADY-STATE EXCITATION 

The use of sinusoidal excitation in the experimental determination of the natural vibra- 
tion characteristics of a system has been a standard technique for many years. The technique 
involves a frequency sweep to determine the approximate locations of the natural frequencies, 
followed by a detailed frequency analysis in the neighborhood of each resonance. During the 
detailed analysis stage, the system is allowed to achieve steady-state vibration at each fre- 
quency before recording the response and moving to the next excitation frequency. Unfor- 
tunately, the method is very time consuming. The slow frequency sweep method was there- 
fore introduced in an attempt to speed up the process. 

In the slow frequency sweep approach the excitation frequency is slowly scanned through 
the range of interest, and it is assumed that the vibration attains quasi-steady-state conditions. 
Choice of the scanning rate depends on the natural frequencies and damping ratios of the 
system under investigation, and the assumption of quasi-steady-state conditions may not be 
valid in many cases. Several errors arise in practice. The response at resonance will be less 
than the steady-state maximum, and the frequency at which the maximum occurs is shifted 
in the direction in which the excitation frequency is changing. Additional errors are intro- 
duced by the averaging time of the function analyzer used to derive the vector diagrams. Back- 
ground noise errors, similar to those encountered in the steady-state sinusoidal method, will 
also be present. Thus, the method has very limited application in practice. 



5.3 TRANSIENT EXCITATION 


5.3.1 Single Pulse 

The frequency response function H(cu), which is the basis for all the damping measure- 
ments, is the Fourier transform of the impulse response function. Thus, if the system under 
investigation is exposed to a unit impulse, the measured response will provide a direct mea- 
surement of H(w) and hence of the damping. In practice the unit impulse is not attainable 
because sufficient energy cannot be supplied to the system in an infinitesimal time period. 
However, an acceptable approximation can be achieved by the use of a single pulse of short 
duration. The use of such pulses has been discussed by White (ref. 25). 

To illustrate the single-pulse method, consider a single-degree-of-freedom system, with 
steady-state frequency response function H(co). If the excitation is p(t) and the response is 
w(t), then 

= < 80 > 

where P(a>) and W(co) are the Fourier transforms of the transient excitation p(t) and response 
w(t), respectively, i.e., 

PM =27jT>(t)e- iwt dt; 

W(co) = ^jHw(t)e‘ lcot dt. 

White (ref. 25) discusses several impulse shapes-rectangular, triangular, and trapezoidal- 
but, for current demonstration purposes, only the rectangular pulse need be considered. Then 

p(t) = A 0<t< T 

= 0 elsewhere, 

where T is the duration of the pulse, and 


P(w) = AT 


sin (<uT/2) -i<oT/2 
(coT/2) 


In certain cases, when the duration of the pulse is short with respect to the shortest 
period present in the vibration of the system (ref. 25), only the response of the system need 
be analyzed without creating large errors in the measured values of the natural frequencies and 
damping ratios. However, it is difficult to determine the magnitude of the errors in advance, 
and in practice usually both the excitation and the response have to be measured. A plot of 
WM/P(co) will represent the frequency response function (equation (80)), and estimates 
of the natural frequency and damping factor can be obtained using the Kennedy and Pancu 
method (ref. 4). 

Although simple in concept, the single-pulse method suffers from several disadvantages 
which limit its usefulness in practice. If the pulse is short, to simulate an ideal impulse and 



excite a wide-frequency range, the energy of excitation and vibration at a particular frequency 
will be low, and accurate measurements will be difficult. Secondly, it is difficult to control 
the frequency content of the excitation, with the result that many modes which are outside 
the range of interest may be excited while modes of interest may not be excited because of 
zeros in the excitation spectrum (see the above example). Thirdly, difficulties may arise in 
maintaining the pulse shape, thereby requiring measurement of the excitation signal as well 
as the response. Finally, noise problems will severely interfere with the accuracy of the results. 


5.3.2 Pulse Train 

As a modification of the single-pulse method, Reed et al. (ref. 23) have investigated the 
use of a series of rectangular and triangular pulses having a continuously decreasing time 
period. Reed et al. claim that the rectangular pulse train has a practical advantage over the 
triangular pulse train and the swept sine wave in that the rectangular pulses are simpler to 
generate. However, the rectangular pulse train suffers from the disadvantage that the energy 
is spread over a large frequency range. Also, pulse trains in general present many of the 
problems associated with the single pulse. 


5.3.3 Rapid Frequency Sweep 

To overcome the difficulties of frequency control, energy input, and pulse shape asso- 
ciated with pulse excitation, an alternative form of transient excitation, the rapid frequency 
sweep, has been introduced by Reed et al. (ref. 23) and Skingle (ref. 24), with further 
developments by White (ref. 14, 25, 26, and 27) and Kandianis (refs. 16 and 28). The rapid 
frequency sweep is achieved by scanning a sine wave signal from frequency cuj to 0 J 2 * n a 
very short time period T. In contrast to the slow sweep, quasi-steady-state conditions are not 
assumed, the response being considered as transient in nature. Obviously there are an infinite 
number of ways by which the frequency can change from cjj to but, for simplicity, a 
linear variation is considered here. Thus, the frequency of excitation can be represented by 
the equation 


to = at + cuj 0 < t < T 
= 0 elsewhere. 

The excitation signal is then 


p(t) = p Q sin (at + wj)dt 

= p Q sin /-^ t“ + cj| tj 0 < t < T 


= 0 


elsewhere. 


(81) 


(82) 


As is the case for the pulse methods, the excitation and response signals are deterministic, 
and the steady-state frequency response function can be obtained by means of equation (80). 



The Fourier transform P(c J) of p(t) is a complicated function, and the characteristics of 
the spectral amplitudes have been summarized by White (ref. 27). The spectrum amplitude is 
not constant within the frequency range toj to 0^2 but h as tw o peaks, one at a frequency of 
cu j + 1.2 *Ja 72 and the other at a frequency ofo >2 - 1 .2 slot/ 2. The mean amplitude of the 
spectrum is7r/2aand the height of each peak is 1.4 times the mean amplitude. Between the 
two peaks, the spectrum is oscillatory about the mean value, the amplitude of the oscillation 
being proportional to \/ yfY. For given values of coj and a long duration time T can be 
selected so that the spectrum amplitude is large enough to avoid problems of signal-to-noise 
ratio and to reduce the spectral ripple to an acceptable value. The cutoff rate of the spectrum 
at a>j and CJ 2 is high. 

Instrumentation to produce the rapid frequency sweep has been developed by White 
(ref. 27), and the modulus of the Fourier spectrum for the excitation shows characteristics 
similar to those described above. The phase spectrum is very irregular. 


In practical applications, the transient method involves problems similar to those 
encountered in the steady-state methods. This is not surprising, since the problems are all 
related to Fourier transformation of the impulse and frequency response functions h (r) and 
H (cS) and to the presence of noise in the excitation and response signals. 

The resolution of close natural frequencies is one example of a common problem. In 
practice this has been interpreted as a requirement on frequency resolution in the power 
spectral density analysis, and as a specification of the maximum time delay for the autocor- 
relation function. For the latter case it was shown (equation (44)) that the displacement auto- 
correlation function for two signals A a sin + 6] and Aq sin [(co^ +co')r + 0] was 


R w( 2 S>t) = (A 2 + A 2 + 2A a Ap cos co'r) 1 / 2 sin (oo d T + 6 +<P(t )) 


for t >0,and 


tan 4>(t ) = 


A d sin co'r 
A a + Aq cos co'r 


(83) 


White (ref. 26) has shown that for the corresponding impulse case, where the total impulse 
function is 


and 


then 


h(t) = hj(t) ±h 2 (t) 
hj(t) = aj sin coit 
h 2 (t) = a 2 sin (go j + co')t, 

h(t) =(a 2 + a 2 +2aja 2 cos co't)*/ 2 sin (cojt + #j(t)), 


(84) 



where 


a 2 sin co't 

tan 0, (t) = — rr ' 

1 7 aj +a 2 cos ou t 

The similarity between equations (83) and (84) is immediately apparent. Positive and 
negative signs in equation (84) are required to take into account the initial phase difference 
(0 or 7 r) between h j(t) and h 2 (t), but this is not necessary in the random excitation case. 

For equations (83) and (84) the maximum time delay has to be sufficiently long to 
identify the frequency w\ This was demonstrated in figures 6 and 7 for the autocorrelation 
function. Under ideal conditions the choice of maximum time delay r m is solely a matter of 
convenience, but in practice it may be dictated by the dynamic range of the analysis equip- 
ment (ref. 26) or by computer storage problems. When severe truncation occurs, the errors 
involved in the estimation of the damping from the measured locus of W(gj)/P(cu) are similar 
to those in the cross-power spectral density (sec. 4.3.5) and single-sided Fourier transform 
(sec. 4.4.4) methods. Correction curves derived from equation (68), examples of which are 
contained in figure 8, can be used to improve the accuracy of the damping estimates from the 
single-sided Fourier transform. 

One important problem, which is difficult to solve by means of the deterministic 
approach to the transient excitation method, is that of noise interference in the excitation 
and response signals. Examples of the distortion in the frequency response function, 
obtained by means of equation (80) when noise is present, are shown in reference 16. Since 
the noise is usually random in nature, Kandianis (ref. 16) has applied random analysis 
methods to the system which, with the exception of the noise, is deterministic in nature. 
Kandianis discusses two possible approaches. The first one computes the power spectral 
density of the response, including the effects of the noise signals. Instead of equation (80) 
for the Fourier components, the power spectral density for the displacement is given by an 
equation similar to that obtained by combining equations (22) and (23). Thus, the displace- 
ment power spectral density function for mode a of the panel is 

S w (x,c4 a = <g(x) |H a M| 2 1 0 2 (x") S p2 (cu) +jA 2 S p| (co) J aa (u) j + S n (w), (85) 

where Sp^Cco) is the power spectral density for the transient excitation applied at x w , Spj(cj) 
is the power spectral density for the noise in the excitation (in the present context this is the 
turbulent boundary layer pressure field), and S n (<x>) is the noise in the response signal. It is 
assumed that Pj(x" t), P 2 (x",t), and n(t) are statistically independent. 

One word of caution is appropriate at this stage. The analysis of Kandianis assumes that 
the excitation signal is stationary and, at a later stage in the single-sided Fourier transforma- 
tion analysis, that the excitation is white noise. This assumption is valid only with certain 
restrictions. Skingle (ref. 24) shows that, under certain assumptions, the autocorrelation 
function for the rapid frequency sweep signal can be written in the approximate form 

sin (co^r) 

r p2M*1^T' < 86) 
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where u >2 is the upper cutoff frequency and the lower frequency cuj is assumed to be zero. 
Equation (86) can be interpreted as implying that there is constant power spectral density in 
the excitation. Referring to the notation of equation (81), the assumptions associated with 
equation (86) are that 


a « 1 and aT^<< 1, 

which in turn imply that » t . Thus, the assumption that the transient signal represents 
white noise is valid only if the maximum time delay in the autocorrelation function is small 
with respect to the duration of the frequency sweep. 

The second approach used by Kandianis (ref. 1 6) is the single-sided Fourier transform 
of the autocorrelation function. The response autocorrelation function obtained from the 
transient excitation in the presence of noise is similar to that encountered in the analysis of 
the autocorrelation function associated with boundary-layer-induced vibration (sec. 4.2.5) 
and suffers from the same experimental problems. Thus, Kandianis proposed the single-sided 
Fourier transform method. The analysis of the method has been discussed in section 4.4 with 
reference to the response to the turbulent boundary layer, and there is no need to repeat the 
analysis in this section. The main difference between the two approaches is that in section 4.4 
the turbulent boundary layer represented the excitation of interest and all other signals were 
noise, whereas in the present case the rapid frequency sweep represents the signal of interest 
and the turbulent boundary layer is classified as noise. 


Combining the results of equations (74) and (75), the single-sided Fourier transform of 
the displacement autocorrelation function will be: 
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The function F w (co) still contains cross terms (a=£ 0) associated with the boundary layer 
spatially distributed excitation. However, the single-point transient excitation does not intro- 
duce additional cross terms, with the consequence that the relative importance of the cross 
terms in equation (87) is less than in equation (74), where only boundary layer excitation is 
present. 



Truncation effects associated with the single-sided Fourier transform, and discussed in 
section 4.4.4, will be applicable to the transient excitation case also. 


5.4 SUMMARY 

The non-steady-state methods discussed in section 5.0 cover a wide range of excitation 
signals. However, the practical value of several of these signals is very limited because of inter- 
ference effects from noise in the excitation or response. All the methods are associated with 
deterministic excitation signals, and there are difficulties in coping with random noise inter- 
ference unless random data analysis methods are adopted. 

For the methods discussed, the rapid frequency sweep method provides the best fre- 
quency control without having to assume quasi-steady-state conditions. In addition, the rapid 
sweep method can be incorporated with random analysis techniques, provided that the maxi- 
mum time delay in the correlation functions is small relative to the time period of the sweep. 
Two steady-state random methods— power spectral density and single-sided Fourier transform 
of the autocorrelation function— have been considered, and the practical problems are very 
similar to those discussed in section 4.0 for steady-state excitation. Within the same constraints 
on time delay, the excitation-response cross-power spectral density method provides a third 
approach for transient excitation, and the problems encountered in practice would again be 
similar to those for the corresponding steady-state excitation method. 



6.0 COMPARATIVE EVALUATION OF METHODS OF DAMPING MEASUREMENT 


In sections 4.0 and 5.0, several methods of measuring the damping of a system have been 
analyzed on an individual basis. Each method is found to have a number of associated prob- 
lems that, in practice, limit the useful range of test conditions. As a consequence, no method 
is applicable to all test conditions. Thus, for a particular experiment, the applicability of each 
method has to be evaluated and the most suitable method(s) selected. An evaluation of this 
sort will be carried out in this section, with specific reference to the test conditions associated 
with the vibration of panels exposed to a moving airstream. 

First it is appropriate to restate the experimental conditions under which the selected 
method of damping measurement will be applied. Damping measurements are to be made on 
rectangular panels that will be mounted in the wall of a wind tunnel. One face of the panel 
will be exposed to turbulent boundary layer pressure fluctuations in the wind tunnel, and the 
other face will be enclosed in a chamber whose volume can be adjusted within a certain range. 
The tunnel will be operated under stable flow conditions. Panel dimensions, at least in the 
initial phases of the experimental program, will not exceed 30.48 cm (12 in.) in length or 
breadth. 

In this section, the damping measurement methods will be compared and evaluated on 
the basis of the theoretical analyses in sections 4.0 and 5.0, with constraints imposed by the 
above experimental conditions. Experimental support for the selected methods will be pre- 
sented in section 7.0, although the test setup in that section will not attempt to reproduce 
the conditions described above. 

Several requirements that the selected methods must satisfy can be identified immedi- 
ately. These requirements are: 

a) The method must be able to exclude, or correct for, background noise in the 
excitation and response signals. 

b) The method must be able to exclude, or correct for, off-resonant vibration in 
neighboring modes. 

c) The method must not require simultaneous measurement of the boundary layer 
excitation and the associated panel vibration, although independent measurements 
may be necessary. 

d) The required instrumentation must be compatible with the space restrictions 
imposed by the test chamber. 

Referring to the analyses in sections 4.0 and 5.0, requirement a) disqualifies the non- 
steady-state methods, with the exception of the rapid frequency sweep method, which uses 
the power spectral density function or the single-sided Fourier transformation of the auto- 
correlation function. The response autocorrelation function for boundary layer excitation is 
eliminated by requirements a) and b), and the excitation-response cross-power spectral density 
method, using turbulent boundary layer excitation, is eliminated by requirement c). A 



combination of requirements a) and d) causes the rejection of the cross-power spectral den- 
sity method when a second excitation is introduced. 

The methods still retained are the power spectral density and single-sided Fourier trans- 
form methods associated with either boundary layer or transient excitation. Retention of the 
power spectral density method may appear surprising at first, since many of the other methods 
have been introduced in the past as improvements over the power spectral density approach. 
Thus, some explanation of the decision is necessary. 

The two main criticisms of the power spectral density method concern the difficulties 
of separating resonant vibration from off-resonant vibration and background noise. Diffi- 
culties that may have been encountered in the past because of too wide a filter bandwidth 
can be adequately overcome with the narrow bandwidths now available through digital analy- 
sis and the use of correction terms. The use of these same narrow filter bandwidths enables 
the separation of close natural frequencies, but there still exists the problem of off-resonant 
vibration and background noise at a natural frequency. Methods can be devised to correct 
for the errors imposed by these interference signals, and one such method has been discussed 
in section 4.0. Although rather crude in form, the correction procedure appears fairly satis- 
factory, particularly if used in conjunction with the optimum positioning of the measuring 
transducer. The more sophisticated methods, which have been proposed as improvements 
over the power spectral density approach, are correspondingly more sensitive to experimen- 
tal scatter. 

The comparison has now reduced to one of steady-state random excitation versus rapid 
frequency sweep. Historically, transient excitation methods were introduced for use in experi- 
mental conditions that are stable for only short periods of time. Since the present test con- 
ditions are stable over relatively long periods, the basic property of the transient excitation 
cannot be classified as an advantage. Further, the method involves the introduction of a sec- 
ond excitation, with accompanying instrumentation problems, and a severe signal-to-noise 
ratio problem because the boundary layer induced vibration is now regarded as noise. The 
transient excitation has to be transmitted to the test panel and, unless a suitable choice of 
panel material has been made, some mass loading of the panel may be necessary. 

However, the transient method possesses advantages associated with single-point exci- 
tation. Mode isolation can be controlled to a certain extent by choice of excitation location, 
although the benefits of the technique are limited by the presence of the boundary layer 
induced vibration. Also, single-point excitation does not introduce cross-term contributions 
in the response spectra, so the adverse effects of the cross terms /3), induced by the 
spatially distributed boundary layer excitation, are reduced. In addition, the rapid frequency 
sweep method uses deterministic excitation, thereby tending to reduce the experimental 
scatter in the measurements. The sinusoidal nature of the excitation allows concentration of 
the excitation energy at selected frequencies and improves the signal-to-noise ratio. 

Basically, the steady-state and rapid frequency sweep methods appear to offer similar 
degrees of accuracy in the damping measurements. Thus, the instrumentation requirements 
will play a major role in determining the most suitable method. Using this as the criter- 
ion, the power spectral density and single- sided Fourier transform methods associated with 
turbulent boundary layer excitation were selected for subsequent experimental demonstra- 



tion. The instrumentation requirements are no different from those normally associated with 
the measurement of structural response to turbulent boundary layer excitation, and the panel 
displacement power spectrum due to boundary layer excitation is an immediate byproduct 
of the damping measurements. The theoretical analysis does not indicate the superiority of 
either the power spectral density method or the autocorrelation single-sided Fourier trans- 
form method, with the consequence that both methods are selected for experimental testing. 



7.0 DATA REDUCTION REQUIREMENTS 


7.1 INTRODUCTION 

Many of the methods of damping measurement discussed in sections 4.0 and 5.0, 
including both of the methods selected in section 6.0, require random process techniques 
for data analysis. It is possible that analog or digital techniques could be used, but, for damp- 
ing measurements, digital methods are recommended because of the greater flexibility in 
choice of filter bandwidth and other characteristics and because the digital methods are 
usually more economical. For some of the Fourier transformation calculations, digital tech- 
niques of data reduction are almost essential. 

The use of digital data reduction methods introduces a number of requirements on data 
sample length, digitization rate, filter bandwidth, etc., and some comments on the require- 
ments are necessary to show their influence on the measured damping. For this reason, the 
present section will provide a general discussion on the requirements. Comprehensive dis- 
cussions of digital data reduction techniques, as applied to random time series, can be found 
in references 20, 29, 30, and 31. 

The data reduction scheme followed in the demonstration experiment, which can be 
considered as a typical procedure, is shown in the flow chart in figure 1 0. The data reduction 
is carried out in four stages that are strongly interdependent, the output of one stage being 
the input of the next. Thus, even the last step in the data reduction cycle can impose restric- 
tions on the data acquisition. As an example, the required statistical reliability of the estimated 
response power spectral density will determine the initial length of the data recording. Although 
some of the data reduction requirements may be peculiar to the current data analysis tech- 
nique, most requirements will be common to all similar damping measurement experiments. 


7.2 DISCRETE FOURIER TRANSFORMS 

A function p(t), continuous in time t, can be represented by a series of discrete samples 
p(n,t) in time. Generally, the discrete points p n are sampled at equal time intervals At, and 
the sampling rate s is defined as s = 1 /At. If N discrete points p n = p(n,Ar) are given for 
n = 1 ,2,...N, then the time series is known in a finite time interval AT where AT = N ■ At, and 
AT can be referred to as a “subinterval.” An infinite discrete time series can be obtained by 
letting N— oo. The complex Fourier amplitude of a continuous function p(t) is given by 
equation (10) as 

p M * ^ jT~P(t)e' iwt dt 

or 

P(f) = f J>(t)e -l27rft dt ( 88 > 
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The equivalent of equation (88) in a finite discrete time series is 

N_i 

, 2 

P k = P(k,Af F )=4 £ P n e- l27rkn / N (89) 

n=-N/2 

where N is an even integer and the fundamental frequency interval 

Af F = ir = N-k7 <90) 

For k = 0, P Q gives the mean of p n . It can be shown that -N/2 < k < N/2, i.e., the 
number of distinct Fourier amplitudes P k is equal to the number N of discrete samples in 
time p n . Thus, AT determines the frequency resolution and s = 1/Ar determines the maximum 
frequency f max = N • Afp/2. The Fourier amplitude P k will have a frequency window of 
the form (sin Afp)/Afp. 

If the direct method of calculating power spectral density is used, as indicated by 
equation (12), then the power spectral density of p n is given by 

S k = S k (k,Afp) = P k P k * (91) 

where the asterisk denotes the complex conjugate. The Fourier amplitude and the complex 
conjugate will each have a frequency window of the form (sin Afp)/Afp. Thus, the spectral 
window of the estimate of the power spectral density S k will be of the form [(sin Afp)/Afp] 2 , 
and has an associated triangular, or Bartlett, lag window. 

When the power spectral density S k (k, Afp) is obtained by means of the direct method 
of equation (91), the autocorrelation function R n = R n (n,Ar) is computed by transforming 
S k (k,Afp). The autocorrelation function will be given in the time interval -N/2 • At < 
t < N • At/2, at discrete time lags At, and the transformation process will incorporate the 
above triangular lag window. 

The single-sided Fourier transform of the discrete autocorrelation coefficients can be 
represented by equation (89). Therefore, the calculated single-sided Fourier transform will 
have a spectral window of the form (sin Afp)/ Afp and the description that follows equation 
(89) will also be applicable. Computer programs that use the fast Fourier transform algorithm 
to calculate power spectral density, autocorrelation, and single-sided Fourier transform func- 
tions will introduce the above characteristics into the computed functions. 


7.3 DATA REDUCTION REQUIREMENTS 

There are several data reduction requirements. Some are peculiar to digital analysis, 
while others are common to all random data analysis techniques— digital or analog. These 
requirements are discussed below. 



7.3.1 Filter Bandwidth 


Errors introduced by finite bandwidths of the analyzing filters are discussed in section 4.0, 
where it is shown that, for rectangular filters, a bandwidth that is half of the measured 
resonance bandwidth, will introduce an error of about 14% in the estimated damping factor. 
Theoretically, corrections can be applied to measured damping factors obtained by means of 
any filter bandwidth, but in practice the correction procedure has acceptable accuracy only 
if the errors are relatively small. Thus, for power spectral density measurements, a filter band- 
width requirement can be specified as Afp < Af^j/2. Obviously, some prior knowledge of 
the damping in the system is required, or an iteration method has to be applied in the 
selection of Afp. 

For a given set of digitized data, the effective filter bandwidth will be a function of the 
time length of the data sample and the digital analysis process used in the computation. The 
frequency-dependent function will be obtained from a time series that can be either the 
original data or a correlation function. Fourier transformation from time to frequency domains 
will be performed with some time or lag window beeause of the finite limits on the time series. 
This lag window will, in turn, determine the characteristics of the spectral window or filter 
in the frequency domain, as indicated in reference 20. For example, a rectangular lag window 
of width AT has a spectral window of equivalent width equal to 1 / AT. Lag windows, such as 
triangular, Hanning, Hamming, and Parzen windows, have the effect of improving the statistical 
reliability of the data but increase the effective bandwidth of the filters for a given value of 
AT. The advantage of these nonrectangular lag windows lies in the reduced amplitude of the 
filter side lobes. 


7.3.2 Total Sample Length 

The length of the data samples used in the random analyses will influence the statistical 
reliability of the reduced data; the shorter the sample length the lower will be the reliability. 
However, long samples require long processing time, and a compromise has to be reached 
between long sample length and practical limitations on cost and computer core storage 
requirements. 

The reliability of the power spectral density estimate is characterized by the normalized 
standard error e where 


1 


(92) 


If the total sample length T consists of V subintervals, each of length AT, then T = V • AT 
and, using equation (90), 



Thus, in practice, statistical scatter can be reduced by increasing the number of sub intervals, 
or subspectra, used to obtain an ensemble average. The effect is illustrated in figure 1 1, where 
the four power spectra are calculated using four different values of V, i.e., V = 9, 25, 50 and 
100 (e= 0.33, 0.2, 0.14, and 0.1). For V=9 (fig. 1 la), resonance peaks in the vicinity of 



1000 Hz are completely obscured by the statistical scatter. As V increases, the statistical 
scatter is reduced and the peaks are better defined. The statistical scatter in the power 
spectral density will introduce equivalent scatter in the autocorrelation and the single-sided 
Fourier transform functions. Thus, when filter bandwidth requirements are specified, based 
on considerations stated in earlier sections, the total sample length requirements can be 
estimated so that a desired degree of statistical reliability can be achieved. 

When the vibration spectrum is obtained by an ensemble average of a series of sub- 
spectra, the vibration autocorrelation function is calculated from the average spectrum by 
means of the inverse Fourier transformation. This method reduces the statistical scatter 
in the autocorrelation function, with a minimum of computer operating time. 


7.3.3 Digitization Rate 

Analog-to-digital conversion of a signal requires the selection of a digitization or 
sampling rate s that is determined by the upper frequency of interest in the spectral analysis. 
Too low a sampling rate will introduce aliasing errors in the spectrum, but too high a sampling 
rate will lead to problems in data handling and economics. Thus, a compromise is again 
required. For Fourier analysis of a time series using the fast Fourier transform algorithm, the 
relation between the upper frequency of interest f max and the lowest acceptable sampling 
rate is given as s = 2f max , but, in practice, a higher sampling rate is necessary. Aliasing errors 
will depend on the spectral content of the signal to be analyzed, but, in general, for stationary 
time series analysis, a sampling rate of s = 2.5f max is adequate. 


7.3.4 Number of Ehta Points 

The filter bandwidth and sampling rate together determine the number of data points to 
be used in the calculations. Thus, the number of data points N is given as N = AT /At = s/Afp 
when Afp = 1/AT. The fast Fourier transform algorithm requires that the number, L, of 
data points used in the transform be an integer power of 2, i.e., L = 1 ) where j is a positive 
integer. If N < L, then (N - L) zeros are added to the data points and the calculation carried 
through with L points. The addition of zeros to the time series will not change the filter band- 
width since this is still a function of the original sample length AT. However, the presence of 
the zeros will change the frequency interval Af' at which the Fourier transform is calculated 
since Af '= 1/AT' where AT' is the sample length after the addition of the zeros. Obviously, 
selection of the number of zeros is again a compromise, with computer storage and cost pro- 
viding a constraint. 

The addition of zeros to reduce the frequency interval has been used in the single-sided 
Fourier transform method (ref. 26) for the analysis of ship vibration. 



8.0 DEMONSTRATION EXPERIMENTS 


8.1 TEST OBJECTIVES 

Two methods— the response power spectral density method and the response autocorre- 
lation single-sided Fourier transform method— have been selected in section 6.0 as being the 
most suitable for the measurement of damping in panels exposed to a moving airflow. These 
two methods have been used to calculate damping factors of a rectangular panel from data 
obtained in laboratory demonstration experiments. The primary objective of the experiments 
is the demonstration of the validity of the two selected methods under laboratory conditions, 
with a secondary objective being the illustration of the various errors and limitations discussed 
in the analytical investigation. Both objectives are essentially met by the demonstration experi- 
ments described in the following sections. 


8.2 DESCRIPTION OF DEMONSTRATION EXPERIMENTS 


8.2.1 Introduction and Planning 

Typical work elements of the demonstration experiments are shown in the flow diagram 
in figure 1 2. The experiments can be divided into three parts: 

a) Fundamental vibration characteristics of the test panel. (Analytical and experimen- 
tal determination of natural frequencies and mode shapes.) 

b) Measurement of panel damping using discrete frequency excitation. (These values 
of the damping are used as standards in the evaluation of the random excitation 
methods.) 

c) Measurement of panel damping using random excitation. 

The test specimen chosen for the experiments was a 30.48- by 17.78- by 0. 127-cm (12- 
by 7- by 0.05-in.) rectangular aluminum panel with clamped boundaries. The 20 lowest 
natural frequencies of the panel were calculated using Warburton’s method (ref. 32), and 
these are listed in table I. It can be seen that the test panel possesses some modes whose 
natural frequencies are fairly well isolated, e.g., modes (1,1), (2,1), and (3,1), and other 
modes with close natural frequencies, e.g., modes (4,2), (5,1), and (1,3) and modes (5,2), 
(3,3), and (6,1). Thus, the two methods of damping measurement can be evaluated for 
single modes and for closely spaced groups of modes. 

As is anticipated from section 3,1, the damping of the test specimen is very low. Since 
this low damping is not necessarily typical of all environments likely to be encountered in the 
wind tunnel, it was considered that the damping measurement methods should be evaluated 
also under conditions of relatively high damping. To this end, the demonstration was carried 
out first on the bare panel and then repeated with a single layer of damping tape (table II) 
applied to one surface of the panel. The damping tape increased the panel damping by an 



order of magnitude and permitted evaluation of the measurement techniques over a wide 
range of values for the damping factor. In the present discussion, therefore, the high and 
low damping cases will refer to the test panel with and without damping tape, respectively. 
The application of the damping has only a small effect on the panel natural frequencies. 

Based on the change in total surface density, the calculated natural frequencies will change 
by a factor of 0.943. 

To make a consistent evaluation of the selected methods, it is essential that there be 
no change in damping of the test panel when the excitation changes from discrete frequency 
to random. This can be achieved with certainty only if the same means of excitation is used 
for both test conditions, and it is for this reason that turbulent airflow was not used as the 
source of the random excitation. There was no guarantee that the acoustic radiation damping 
would not change in the presence of the airflow. 

Two possible excitation methods, single point or distributed loading, were available for 
the experiment. Single-point excitation has the advantage that selected modes can be excited 
or suppressed. However, since the panel material was aluminum, noncontacting magnetic 
excitation could not be used and single-point excitation would have involved mass loading 
of the panel. Acoustic excitation does not have the modal selectivity of single-point loading 
but is more typical of the loading distribution imposed by a moving airstream. Furthermore, 
the two selected measurement methods do not specifically require point excitation. Thus 
acoustic plane wave excitation was chosen for the demonstration experiments. 

Knowledge of the panel mode shapes is highly desirable, even essential, if reasonable 
accuracy is to be achieved in the damping estimates. Identification of the nodal and antinodal 
lines allows selection of the measuring transducer locations so that vibration in a particular 
mode can be given preference. Qualitative information, such as that obtained from Chladni 
or nodal patterns, is sufficient for the selection of the measurement positions. 


8.2.2 Experimental Setup and Instrumentation 

A block diagram of the experimental setup and instrumentation is shown in figure 13. 
The aluminum test panel was mounted on a heavy pedestal such that it could be excited by 
acoustic waves with angles of incidence varying from grazing to normal by changing the 
orientation of the acoustic horn and panel. The test panel could be mounted in the vertical 
or horizontal plane, the horizontal plane being used only when obtaining the Chladni figures. 
The excitation was either single-frequency (pure-tone) or band-limited white-noise acoustic 
waves, and the intensity of the excitation was controlled by means of the voltage across the 
acoustic horn. 

The panel response was measured in terms of displacement, using a noncontacting Photocon 
proximity transducer, Model PT5, which was mounted on a traverse bar. It was possible to 
locate the probe at any position over the plate, with the gap between the probe and the test 
panel adjusted to any desired value. The gap between the probe and the test panel at rest was 
chosen as 0.076 cm (0.03 in.) to provide the required sensitivity for the displacement meas- 
urements. The Photocon proximity probe was calibrated, and the ac voltage output of the 
probe provided a measure of the displacement of the panel. 



The excitation signal could be measured either in terms of the voltage input to the driver 
unit of the acoustic horn or in terms of the acoustic output of the horn. To measure the sound 
field, a 1 /4-inch-diameter Bruel and Kjaer microphone was located near the horn exit. When a 
band-limited white-noise signal was applied to the driver unit, the acoustic pressure field in 
the neighborhood of the test panel had a power spectral density that was flat to within ±5 dB 
in the frequency range of 200 to 2000 Hz and ±2 dB in any octave band in this frequency 
range. The horn response did not contain sharp spikes but varied slowly with frequency. In 
the vicinity of any panel natural frequency, the horn response could be considered to be flat. 

During the discrete frequency tests, the phase angle between the voltage input to the 
horn and the voltage output of the proximity or displacement probe was measured in addition 
to the amplitudes of the two signals. The data were analyzed on line. For the random excita- 
tion methods, the voltage input to the driver, the acoustic output of the horn, and the panel 
response were recorded on magnetic tape to provide an analog source tape. The analog signals 
were subsequently passed through an analog-to-digital converter to produce a digital data tape. 
A CDC 6600 digital computer was used to process the digitized data. 

The experiments were conducted in an acoustical test room where the acoustic environ- 
ment could be controlled. 


8.3 EXPERIMENTAL PROCEDURE 


8.3.1 Preliminary Experiments 

The preliminary experiments included experimental determination of the natural fre- 
quencies, qualitative estimates of mode shapes by means of Chladni figures, determination 
of suitable angles of incidence of the acoustic waves for exciting desired modes, and coarse 
estimates of the half-power bandwidths of a few modes. Using acoustic waves of grazing inci- 
dence, it was possible to excite all mode orders in the direction of propagation of the acoustic 
waves. Therefore, waves of grazing incidence propagating along the length of the panel were 
used in most of the investigation. Only when obtaining the Chladni figures for modes with 
even order along the width of the panel was the excitation changed such that waves of grazing 
incidence would propagate along the width. 

To obtain the nodal line patterns (Chladni figures) for the test structure, the panel was 
mounted horizontally and white, dry sand or fine aluminum filings were spread thinly on the 
panel surface. The panel was excited at each natural frequency in turn, and the nodal pattern 
was photographed when the excitation was turned off. 

Figure 14 shows Chladni figures for the bare panel. Chladni figures for the panel with 
damping tape are almost the same as the corresponding figures for the panel without tape. 
Thus, the patterns shown in figure 14 can be considered as representative of the low and high 
damping conditions. The irregular nodal patterns that occur at some natural frequencies are 
most likely due to inhomogeneities in the thickness and flatness of the panel and to initial 
stresses. However, for the present demonstration, these details are irrelevant. 



One peculiarity observed during the preliminary experimental stage is worthy of com- 
ment. In the process of obtaining a quick estimate of the half-power bandwidth of the (2,1) 
mode for the low damping case, the response amplitude varied significantly even when the 
amplitude and frequency of the excitation were held constant. A systematic investigation of 
the phenomenon ruled out the possibility of drift in the measuring instruments, and it was 
deduced that the variation in the amplitude was inherent to the panel response and was due 
to drift in natural frequency. The frequency drift (measured via variation of response ampli- 
tude as indicated above) showed a strong dependence on the peak response amplitude. When 
the panel displacement was slowly increased, the first noticeable occurrence of the drift was 
observed at a peak amplitude of about 0.0102 mm (0.0004 in.). The drift increased signifi- 
cantly when the peak displacement exceeded 0.025 mm (0.001 in), but no drift was observed 
below 0.0076 mm (0.0003 in.) peak displacement. Therefore, the excitation was adjusted to 
give a peak displacement of less than 0.0076 mm (0.0003 in.). A similar, but less pronounced, 
frequency drift was observed for mode (1,1). No frequency drift was observed for the high 
damping case. 


8.3.2 Discrete Frequency Excitation 

For this phase of the experiment, similar procedures were used for the high and low 
damping cases. The panel was excited by a pure tone at the natural frequency of a mode 
under investigation. The measuring displacement probe was located at an antinode of the 
mode, the location being selected by means of the previously obtained Chladni figures. The 
excitation was then adjusted such that the peak deflection at resonance was less than 0.0076 mm 
(0.0003 in.). For the higher order modes and for the high damping case, the maximum allow- 
able input voltage to the driver was often the limiting factor; under these conditions, the peak 
deflection was substantially less than 0.0076 mm (0.0003 in.). The response amplitude and 
the phase with respect to the excitation were measured at a series of frequencies in the vicinity 
of the resonance. For a few selected modes, the excitation frequency was varied in both 
ascending and descending directions to ensure repeatability. 

In the case of a mode with a natural frequency close to other natural frequencies, the 
contributions from other modes were minimized by locating the measuring probe on or near 
nodes of the undesired modes. In such a case, the probe may be located off the antinode of 
the mode under investigation. For example, when investigating mode (6,1), the measuring 
probe was located at xj = 8.89 cm (3.5 in.), X 3 = 1 1.43 cm (4.5 in.) so that the contribution 
of the (3,3) mode was eliminated, the probe being located on a node line of the (3,3) mode 
(see fig. 14), and the contribution of the (5,2) mode was minimized. Similarly, when investi- 
gating the (5,2) mode, the contribution of the (3,3) mode was eliminated by locating the 
probe at xj = 12.07 cm (4.75 in.), X 3 = 13.34 cm (5.25 in.). Measurement locations used 
during the discrete frequency test are given in table III. 

The effects of background noise and off-resonant contributions have been discussed in 
section 4.0. To obtain estimates of the magnitude of these contributions in the displacement 
spectra, the response spectra were investigated on either side of each natural frequency, and 
the spectral minima adjacent to each natural frequency were identified. The spectral densities 
and frequencies associated with these minima were recorded. 



8.3.3 Random Excitation 


A procedure similar to that described above for discrete frequency excitation was used 
for the random excitation, except for the following two differences: (1) band-limited white 
noise was substituted for the pure tone and (2) the excitation, the response, and the input 
voltage to the driver of the acoustic horn were recorded on a magnetic tape to produce an 
analog source tape. Measurement locations and excitation bandwidths are given in table IV. 

Reduction of the random data for the demonstration experiment was carried out in 
four steps: data acquisition and recording to produce the analog tape, analog-to-digital con- 
version, use of a computer program to calculate the displacement power spectral density and 
autocorrelation functions, and, finally, calculation of the single-sided Fourier transform of 
the autocorrelation function. The procedure is the same as that outlined in figure 10. 

The voltage signals representing the panel displacement, acoustic excitation, and current 
supplied to the driver unit of the acoustic horn were recorded on magnetic tape at a tape 
speed of 152 cm/sec (60 in. /sec), to produce the analog source tape. The analog tape was 
then replayed into digitizing equipment to produce a digital tape. Since the maximum fre- 
quency of interest was specified as 2000 Hz, a digitization rate of 8000 samples per second 
was used. This satisfied the condition s >4f max , required by the single-sided Fourier trans- 
form computations (see below), for all run conditions. Digitizing rates less than 8000 samples 
per second could be used if f max < 2000 Hz, and these lower sampling rates were obtained 
from the original digitized data by means of a decimation process. However, the condition 
s < 4f max was always satisfied. A total of 60 sec of data was recorded for each test run and, 
from this, a total of 54 sec was digitized. 

For the high damping case, all data runs were averaged over 50 subspectra. Thus, the 
normalized standard error e of the power spectral density estimates was about 0. 1 5 for the 
high damping case. When the low damping data reduction was carried out, the maximum 
available sample time length of 54 sec became a critical factor, and in only three runs (7A, 
7B,and 8B) could the spectra be computed by averaging 50 subspectra. A smaller number of 
subspectra was used in all other low damping runs. 

Filter bandwidth requirements are extremely critical at low frequencies in the low 
damping case. For example, in modes (2,1) and (3,1), filter bandwidths of 0.35 and 0.9 Hz 
are desirable. However, the limited data sample length of 54 sec was not sufficient to provide 
acceptable statistical scatter (e < 0. 1 5) for such narrow filter bandwidths, so a larger band- 
width of about 1.0 Hz was used. Filter bandwidth did not pose a problem for the high damp- 
ing results, the condition Afp < Af^j/2 being satisfied in all cases. 

The computer program used to calculate the single-sided Fourier transform imposed 
additional constraints on the data. The program was checked with respect to known analyti- 
cal functions and it was found that, for a transformation error of less than 1%, Afp should be 
less than Afj^j/4 and digitizing rate s should be greater than 4f max . For an error of less than 
15%, Afp < Afjyj /2 and s > 4f max . In the computation of the single-sided Fourier trans- 
form, zero terms were added to the autocorrelation function so that an adequate number of 
data points was available for the plotting of the response function. Addition of zero terms 
increased the maximum time delay by a factor of four for the low damping case, and by a 
factor of at least seven for the high damping case. 



8.4 DISCUSSION OF RESULTS: DISCRETE FREQUENCY EXCITATION 


Three possible ways of calculating damping factors using the discrete frequency excita- 
tion are presented in figures 15a, b, and c. For a single mode, such as the (2,1) mode, which 
has a natural frequency well separated from the adjacent natural frequencies, all three meth- 
ods would yield damping estimates of similar accuracy. Because of the low damping in the 
case of the bare test panel, most of the modes can be considered to satisfy, approximately, 
the conditions for a single mode. Thus, for the low damping case, the response amplitude 
method and the Kennedy-Pancu vector diagram (amplitude phase plot) give damping factors 
that are within an acceptable experimental scatter, e.g., the damping factors calculated using 
the two methods are within ±2.5% of each other for all modes except three. Two of the three 
exceptions have experimental scatter of ±5%, whereas, for the remaining mode (5,2), the 
Kennedy-Pancu method gives a damping factor that is 67% of the one calculated from the 
response amplitude method. The higher estimate from the response amplitude method is 
attributed to the off-resonant contribution from the (6,1) mode. 

When damping tape is applied to the test panel, only a few lower order modes, e.g. (2,1), 
(3,1), and (4.1), satisfy the conditions for a single mode. Due to contributions from the off- 
resonant modes, damping factors calculated from the response amplitude method are gener- 
ally higher (by an average of 15%) than those calculated by the Kennedy-Pancu method. No 
background noise corrections are applied to the discrete frequency response amplitude data 
because of problems in determining the correct phase-angle relationships. 

The effect of background noise and contributions from the off-resonant modes can be 
readily seen from the Kennedy-Pancu diagram. For a single-degree-of-freedom system, when 
the off-resonant contribution is zero, the center of the resonance circle will be located at half 
the peak amplitude, i.e., at 0.5 when the normalized amplitude is used. Thus, a deviation from 
0.5 is a measure of the off-resonant contributions. As an example, figures 16 and 17 show the 
Kennedy-Pancu diagrams from the low and high damping cases, respectively, of mode (2,3). 
The centers of the resonance circles for the low and high damping cases are located at 0.44 
and 0.26, respectively. 

As a contrast, consider mode (2,1), which is well separated from other modes. The 
Kennedy-Pancu diagram for the (2,1) mode, even for the high damping case (fig. 18), shows 
that the center of the resonance circle is at 0.46. 

When the natural frequencies of the two modes are so close that the two resonance 
circles cannot be clearly defined in the frequency range of the half-power bandwidth of the 
response peak, the damping factors can be calculated from the phase-angle plot (fig. 15b). 

The damping factors for the (5,1), (4,2), and (1,3) modes for the high damping case were 
obtained using this technique. All other damping factors shown in table III were obtained 
using Kennedy-Pancu diagrams. The damping factors shown in table III are used as the datum 
damping factors for the test panel with and without damping tape. 



8.5 DISCUSSION OF RESULTS: RANDOM EXCITATION 


The stochastic equivalents of the three schemes of calculating damping factors are shown 
in figures 1 5d, e, and f. Once again, in principle, all three methods can yield estimates of sim- 
ilar accuracy for a single mode that is free from background noise and contributions from off- 
resonant modes. The practical application of these methods to a multimodal system is 
described here. The primary objective is to assess the response power spectral density and 
single-sided Fourier transform of the response autocorrelation methods. Therefore, no 
attempt is made to calculate the damping factors of each and every mode. Since the panel 
was excited by acoustic waves propagating along the length of the panel, most of the modes 
with mode order higher than one in the width will not be excited very efficiently. Those 
modes are usually excluded from the discussion, but this exclusion does not cause any loss 
of generality in the conclusions. 


8.5.1 Response Power Spectral Density Method 

The high damping case, which will be affected by the background noise and contribu- 
tions from off-resonant modes more than the low damping case, is considered first. The 
response spectra calculated from runs 2B, 4B, and 9B (see table IV) are shown in figures 19, 

20, and 21, respectively. Spectra drawn in figures 19 and 20 indicate that the peak response 
in the (2,1) and (4,1) modes is higher than the nearest valley in the spectrum by a factor of 
more than 10. Thus, these two modes are considered to represent a single-mode case. For 
both the modes, Af^/Afp >5; therefore, errors due to filter bandwidth are less than 2%. 

The calculated damping factors for the (2,1) and (4,1) modes are 1.08 and 1.12 times the 
datum damping factors (see table V). 

Now consider mode (3,1) shown in figure 20. Due to the measurement location for run 
4B (see table IV and Chladni figures in fig. 14a) the (3,1) mode is suppressed and the response 
peak shows some contribution from the background noise. The damping factor for mode 
(3,1), before applying background noise correction, is 1 .28 times the datum damping factor. 
When the background noise correction is applied (fig. 4), the corrected damping factor in 
table V for the (3,1) mode is 1.11 times the datum value. 

A relatively large excitation bandwidth, 200-2000 Hz, was used for run 9B. Thus, the 
response spectrum shown in figure 2 1 will be similar to that obtained by exciting a panel by 
broadband turbulent boundary layer pressure fluctuations. The measurement location was 
such that the probe would measure all modes with mode order less than six in the length 
direction and less than four in width direction. As expected, the spectrum shows a series of 
peaks with different relative levels, and the (3,2) mode can be taken as an example of very 
high contributions from background noise. For this mode, the damping factors calculated 
with and without background noise correction are, respectively, 1.06 and 1.33 times the datum 
damping factor. 

Measured damping factors, estimated by means of the power spectral density method, 
are shown in table V for several modes in the high damping case. Where necessary, the meas- 
urements have been corrected for background noise using figure 4. The data are separated 
into two categories that are determined by the magnitude of the background noise. For 



table Va, the displacement power spectral density ratio of resonant peak to adjacent valley 
is greater than 10 dB, whereas, in table Vb, the ratio is less than 10 dB. In all cases, the filter 
bandwidth is less than a quarter of the measured resonance bandwidth so that the resolution 
error, predicted by figure 5, is less than 3%. For the high damping case, when the background 
noise is low, table Va shows that the error in the measured damping factor varies from -11% 
to +25% of the corresponding datum value, with a root mean square (rms) error of 12%. 
When the background noise is high (table Vb), the error increases and, for the data shown, 
lies in the range of +6% to +89% with an rms error of 44%. 

The problems encountered in the measurement of low damping factors are rather dif- 
ferent from the high damping case. The off resonant/background noise contributions do not 
now pose problems. However, due to the narrowness of the resonance peaks, the filter band- 
width becomes a critical factor for the lower order modes. Also, as discussed in section 7.0, 
long sample lengths are required for good statistical reliability with narrow filter bandwidths. 
The response spectra calculated from low damping runs 4B, 9B, and 8B are shown in figures 
22, 23, and 24, respectively, the total sample length being the same for all three cases. Filter 
bandwidth Afp, is 0.5 Hz in figures 22 and 23 and 1.0 Hz in figure 24. The two cases with 
the smaller filter bandwidth show higher statistical scatter, which makes it difficult to calcu- 
late accurately the half-power bandwidths of the peaks. 

A comparison of the spectra in figures 23 and 24 indicates another important problem 
area. The excitation bandwidths for runs 8B and 9B are 800-1600 and 200-2000 Hz, respec- 
tively. Due to the wider excitation bandwidth for run 9B, the lower order modes, which can 
be excited relatively efficiently, dominate the spectrum in figure 23. In run 8B, the excita- 
tion has an 1 8-dB per octave rolloff below 800 Hz, thus, the lower order modes are de-empha- 
sized in figure 24, and the peak response amplitudes are of similar magnitude for many of the 
modes shown. When the panel is excited by the broadband noise due to turbulent airflow, 
an effect similar to that in figure 24 can be obtained by filtering the response signal from the 
displacement probe before recording the data. 

Measured damping factors for the low damping case are shown in table Vc for the power 
spectral density method. In all cases except one, the power spectral density peak-to-valley 
ratio is greater than 10 dB, and the damping factor lies within a range of -18% to +32% of the 
datum values. The rms error is 15%. Except where indicated, the damping factors in table Vc 
were obtained using filter bandwidths less than half of the measured peak bandwidth. 


8.5.2 Single-Sided Fourier Transform Method 

The single-sided Fourier transform F w (to) of the displacement autocorrelation function 
was computed using the fast Fourier transform algorithm on a digital computer. Data reduc- 
tion requirements associated with the computer program have been discussed in section 8.3.3. 
These requirements were satisfied for all modes investigated at the higher damping condition 
but not for all modes with the lower damping. 

Consider first the high damping case because these results show the accuracy of the 
method when used under satisfactory conditions. As an example, the function F w (cj) for the 
(2,1) and (2,3) modes is shown in figures 25 and 26, respectively. In both cases, the locus of 



F w (c o) is approximately circular with some distortion because of truncation errors introduced 
by the finite maximum time delay. The curves can be compared with the corresponding fig- 
ures for discrete frequency excitation shown in figures 18 and 17, respectively. Vector curves 
for modes (5,1), (4,2), and (1,3), which have close natural frequencies, are shown in figure 
27. The figure clearly shows the effect of contributions from off-resonant vibration. 

In common with the function H (co), the natural frequency of a mode is located where 
the rate of change of arc length with frequency is a maximum for the function F w (e o). In 
many cases, such as figures 25 and 26, the natural frequency can be located by visual inspec- 
tion of the vector diagram. The resonance diameter is drawn and the damping factor calcu- 
lated as illustrated in figure 15f. When there are close natural frequencies, it may be difficult 
to locate each resonance by visual inspection of the vector diagram, and a curve of arc length 
as a function of frequency is plotted as shown in the example of Figure 28. If the resonance 
circle is not clearly defined over the half-power bandwidth, the damping can be calculated 
from the phase-angle plot in the neighborhood of the natural frequency. 

Measured damping factors for several modes, estimated using the single-sided Fourier 
transform of the displacement autocorrelation function, are shown in table VI. Data for the 
high damping case have been separated into low and high background noise categories, as in 
table V for the power spectral density method. For modes with small background noise and 
off-resonant contributions, the estimated damping factors in table Via are within -18% to 
+23% of the datum values, with an rms deviation of 13%. For high-noise conditions, defined 
as a power spectral density peak-to-valley ratio less than 10 dB, the estimated damping factors 
show deviations within the range -1 1% to 77% of the datum values with an rms error of 47%. 

One of the requirements of the computer program for calculating the single-sided 
Fourier transform is that the frequency resolution of the transform would be less than a 
quarter of the half-power bandwidth of the resonance peak (sec. 8.3.3). Within the limits of 
the data reduction system used in the demonstration experiments, the above requirement 
could not be satisfied for all lightly damped modes, such as modes (2,1) and (3,1) in table Vic. 
The damping factors in table Vic, which refer to the lightly damped conditions, have been 
corrected for truncation errors based on figure 9 but still show a wider experimental scatter 
than for the higher damping. For the data in table Vic, the estimated damping factors differ 
from the datum values by amounts ranging from -4% to +70%, with an rms deviation of 43%. 


8.5.3 Comparison of Random Analysis Methods 

Damping factors for several normal modes of a rectangular panel have been estimated, 
from measurements made under random excitation conditions, using displacement power 
spectral density and single-sided Fourier transform methods. A comparison of the results 
identifies two important parameters, one of which is the presence of background noise and 
the other is the frequency resolution requirement of the single-sided Fourier transform method. 

Under the best test conditions, with low background noise such that the displacement 
power spectral density peak-to-valley ratio exceeded 10 dB, and with good frequency reso- 
lution defined as being a filter bandwidth less than a quarter of the resonant peak bandwidth, 
the measured damping factors were close to the datum values. The rms errors were about 13% 



for both of the methods tested. When the background noise or off-resonant vibration 
increased but the high-frequency resolution was maintained, the measured damping factors 
were less reliable, but both methods were of the same accuracy (an rms error of about 45%). 

Differences between the two methods were observed when the resolution criterion was 
not satisfied, even though an adequate signal'-to-noise ratio was achieved. This is shown in the 
damping estimates for the low damping data. The power spectral density method shows an 
rms error of about 1 5%, which is close to the 1 3% error in the high damping data, but the single- 
sided Fourier transform method has an rms error of about 43%. It is apparent that the fre- 
quency resolution requirements imposed by the single-sided Fourier transform computation 
are more stringent than those for the power spectral density method. In certain test cases, 
with high damping and close natural frequencies, no damping estimate could be obtained 
from the power spectral density function, but an estimate could be obtained using the single- 
sided Fourier transform. However, this estimate could be highly inaccurate. 


9.0 CONCLUSIONS AND RECOMMENDATIONS 


Several methods of measuring the damping of panels have been evaluated with particular 
emphasis being placed on the practical problems associated with the use of the methods. In 
particular, the techniques were evaluated for application to the measurement of the damping 
of panels exposed to turbulent airflow. 

Based on analytical analyses, several methods were judged to be unsuitable for the 
proposed application. Non-steady-state methods, which do not incorporate random data 
analysis techniques, were rejected because of noise problems, and the response autocorrela- 
tion method was rejected because of problems of noise and off-resonant vibration interference. 
The excitation-response cross-power spectral density function, using either boundary layer or 
a second excitation, was eliminated either because it required simultaneous measurement of 
the boundary layer pressure field and the vibration, or because of potential low signal-to-noise 
ratios and problems of instrumentation location in the pressure equalization cavity. 

Two basic methods were retained: the response power spectral density and the single- 
sided Fourier transform of the response autocorrelation function. For these methods, two 
alternative excitations, the turbulent boundary layer or a rapid frequency sweep, were con- 
sidered. However, the latter excitation was rejected because of additional instrumentation 
requirements and because there was no need for transient techniques. 

The power spectral density and single-sided Fourier transform methods were demon- 
strated experimentally under conditions of spatially distributed, random excitation. Under 
many test conditions, the two methods estimated damping factors with similar experimental 
accuracy. However, the single-sided Fourier transform method imposes more stringent require- 
ments on frequency resolution than does the spectral density method. This becomes apparent 
when the damping is very low and computer restrictions do not permit the resolution require- 
ment to be satisfied. On the other hand, the single-sided Fourier transform method has advan- 
tages when the damping is high and natural frequencies are close together. Under such con- 
ditions, the power spectral density method has little value. 

From this evaluation, it is recommended that, in general, the response power spectral 
density method be used to measure the damping of panels exposed to turbulent airflow. Pro- 
vided appropriate correction terms are introduced, the reliability and accuracy of the method 
appears to be as good as that associated with other methods, and the data acquisition and 
reduction requirements are no more than those normally demanded for structural response 
measurements. However, if the particular conditions of high damping and close natural fre- 
quencies are important, then the single-sided Fourier transform of the response autocorrela- 
tion function should be used for improved accuracy. 


Commercial Airplane Group 
The Boeing Company 

Seattle, Washington, August 1971 



APPENDIX A 


Finite Bandwidth Effects in the Autocorrelation Method 


The autocorrelation function for a band-limited response signal can be studied by means 
of the function 


R w (t;Aco) = J[ <Mco)|H a (w)| 2 e lcJT dw, 


where 


< J >( to )= 1 w i < M < w 2 

= 0 elsewhere 

Using convolution relationships, equation (Al) can be written in the form 


where 


R w (t;Auj) = ^rjf~H'(tWr ' t)dt 


H'(t) = j[°°|H a (co)| 2 e icot dw 


(Al) 


(A2) 


(A3) 


(A4) 


and 


0(t)= f°°< i>(co)e itot dw. 

•'-oo 


From equations (A2) and(A5), 
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<t>{ t) = -^(sinc^t - sincj^t), 


and from equation (A4) and definition of H a (iu), 
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The parameter is the frequency of free damped vibration, $Ct ) . Substitut- 

ing equations (A6) and (A7) in (A3), 


R w (t;Aco) = R w (t;A«) w - R w (t;Aw) w ] 


where 


R w (t;Aco) w 


2 2M 


1 /*<*> -OOvf/vt / 

2 3y X e cosw d t + — ^smaidt 

a w a ^ \ (1 - f a 2 ) 


sincu 7 (r - t) sincoofr + t)\ 
(T-t) + (r + 1) ) dt 


Two approaches are possible to solve equation (A8). In one approach, it can be shown 


t f°°-at i* sinc(r+t), sinc(T-t) J4 
1, = J e at c°sbt I ( / +1) + (f . t) jdt 


/ c coszr j_ 
-5 =-dz 

-c a 2 + (b + z) 2 


>2-r 


g-at sinht sin c(t + t) sinc(r-t) 
e sinot (T+t) i- (T _ t) 




c z coszr 
c a 2 + (b + z) 2 


The integrals in equations (A9) and (A10) can be related to tabulated exponential 
integrals for numerical evaluation. However, a better insight may be possible using the 
approach of Kandianis (ref. 15). Equation (A8) contains integrals of the form 

sinco 2 (r±t) 

X e coscudt (r±t — dt 


r°°-H*r a t. sinco 2 (r±t) 

X e smw d t_ (r±“t) dt 
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As an example, select 



e cosco^t 


sinciWr + 1) 

(r+t) 


dt 


^afa^f 00 ~W/vf a t _ sinc^t r T -W/Y? a t , sina^t 
= e 01 J e a 01 cosu)j(t - r) — j-=-dt - J e a costu^t - t) — - — dt 


(All) 


The infinite integral on the right side of equation (All) can be evaluated using the results of 
Kandianis (ref. 15), 


if(w a 2 - w 2 2 )>0 

+ n if (to a 2 - 0^2^) < 0 
and 

(r a QJ a ) 2 + Ht + t0 2^ 2 
%“a )2 + Hx ’ w 2> 2 _ 

— CO t | T | 

Terms of this type provide constant factors of the exponential e 01 a and so do not 
influence the damping estimate. 

To evaluate the finite integral on the right hand side of equation (All), Kandianis (ref. 
15) makes use of the reduction formula 


if 


°°-War a t. sinco 2 t t 
e sinw d t — — dt = -log e 


e cosco d t 


— r~dt =4tan- J 



J*xPe ax sin bxdx = -77-^ — ^(asinbx + bcosbx ) — ^ — x fxP ^e ax (asinbx - bcosbx)dx (A12) 
a 2 + b 2 a 2 + b 2 

Taking the first term on the right side of equation (A 1 2) as being an approximation for the 
integral on the left side, Kandianis shows that the autocorrelation function for a band-limited 
signal contains terms of the form (cos a>2T)/rand (sin o>2 r)/r, with similar terms in cjj. These 
terms, since they are functions of r, will distort the autocorrelation function and introduce errors 
in the damping estimates. Thus, the filter bandwidths have to be chosen to minimize these 


errors. 



APPENDIX B 


Truncation Effects in the Single-Sided Fourier Transform Method 


In practice, truncation errors arise in the computation of the single-sided Fourier trans- 
form function, and it is necessary to correct the measured damping factor. The magnitude of 
the correction can be estimated by means of the integral 


where 


F w<“;%> - r, j[“D(r)R w (r)e-i“'dr, 


D(t) =1 0 < r < t 


m 


= 0 t>t 


m 


Using the faltung relationship, equation (Bl) can be written in the form 


(Bl) 
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S W M =^S_Z R ^ T)e ~ i0)TdT 

= S p (w)|H a (w)| 2 


(B4) 
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for a single-degree-of-freedom system with frequency response function H a (to). Substituting 
(B4) and (B5) in (B3) gives 
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where it is assumed that the excitation is white noise so that Sp(co) = Sp. Substituting 
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in equation (B6), and using contour integration in the upper half plane, equation (B6) becomes 


Z 7 ^ y 2 t x ( w ) + iyMl 

4M a 2 co a 4 f a (l - f a 2 ) 
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If (l/p)(ds/dw) is taken as a convenient measure of the characteristics of the locus of 
F w (co;T m ), as suggested by Clarkson and Mercer (ref. 1 7 ), then 
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Evaluating (l/p)(ds/dco) at the frequency of maximum value of |H a (co)| 2 , i.e., at 

w r 2 = d - 2 r a 2 )^ a 2 

and assuming that the damping is small so that f a 2 « 1 , equation (B8) reduces to 

d 2 x 
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or, after substituting for d^x/dco^ and dy/dco, 


_[ds _ 47T 
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TABLE I. -CALCULATED 9 NATURAL FREQUENCIES OF 
ALUMINUM PANEL b WITH CLAMPED EDGES 


m c 

1 

2 

3 

4 


Frequency , d Hz 

1 

253 

642 

1226 

2004 

2 

359 

742 

1327 

2106 

3 

541 

913 

1492 

— 

4 

798 

1157 

1728 

— 

5 

1124 

1473 

2035 

— 

6 

1519 

1861 

— 

— 

7 

1980 

— 

— 

— 

8 

— 

— 

— 

— 


a Warburton method (ref. 32) 

b 30.48x 17.78 x 0.127cm (12 x 7 x 0.05 in.) 

c Mode order (m,n) defined by number of half wavelengths along 
length and breadth of panel. 

^Dashes indicate natural frequencies higher than 2150 Hz. 


TABLE II.— DETAILS OF DAMPING TAPE 


Type 

Scotch brand 428A 

Thickness of aluminum foil 

14.0 x 10' 3 cm (5.5 x 10' 3 in.) 

Thickness of adhesive 

6.4 x 10' 3 cm (2.5 x 10' 3 in.) 

Thickness of foil plus adhesive 

20.3 x 10' 3 cm ±10% (8 x 10' 3 in. ±10%) 

Weight 

4.40 N/m 2 (0.09 Ib/sq ft) 

Weight of bare panel 
Weight of panel and tape 

0.89 
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TABLE I II. -PAN EL NATURAL FREQUENCIES AND 


DAMPING factors measured with discrete frequence excitation 


Mode 

order, 

Location of Proximity Probe 

High damping 

Low damping 

X 

1 


x 3 

f _ ,Hz 
m,n 

^*m,n 

f n ,Hz 
m,n' 

5 m,n 

m,n 

cm 

in. 

cm 

in. 

2,1 

7.62 

3.0 

8.89 

3.5 

359 

0.0153 

362 

0.00102 

3,1 

15.24 

6.0 

8.89 

3.5 

538 

0.0115 

545 

0.00170 

1,2 

15.24 

6.0 

13.34 

5.25 

599 

0.0048 

612 

0.00202 

2,2 

7.62 

3.0 

13.34 

5.25 

700 

0.0089 

717 

0.00163 

4,1 

11.43 

4.5 

8.89 

3.5 

773 

0.0175 

778 

0.00190 

3,2 

15.24 

6.0 

13.34 

5.25 

867 

0.0113 

888 

0.00127 

5,1 

15.24 

6.0 

8.89 

3.5 

973 

0.0200 

1095 

0.00127 

4,2 

15.24 

6.0 

8.89 

3.5 

1057 

0.0163 

— 

— 

1,3 

15.24 

6.0 

8.89 

3.5 

1098 

0.0127 

1124 

0.00188 

2,3 

7.62 

3.0 

8.89 

3.5 

1184 

0.0127 

1221 

0.00313 

3,3 

15.24 

6.0 

8.89 

3.5 

1351 

0.0097 

1405 

0.00194 

6,1 

8.89 1 

3.5 

11.43 

4.5 

1351 

0.0138 

1404 

0.00164 

5,2 

12.07 

4.75 

13.34 ' 

5.25 

1355 

0.0137 

1405 

0.00231 


TABLE 1 1 /.-MEASUREMENT LOCA TIONS AND 
EXCITATION BANDW/DTHS FOR RANDOM EXCITATION 


Run 

Location of proximity probe 

Frequency range, 3 

Hz 

x i 

x 3 

cm 

in. 

cm 

in. 

2B 

7.62 

3.0 

8.89 

3.5 

200-400 

3A 

15.24 

6.0 

8.89 

3.5 

400-800 

3B 

7.62 

3.0 

13.34 

5.25 

400-800 

4A 

5.08 

2.0 

13.34 

5.25 

200-800 

4B 

19.05 

7.5 

13.34 

5.25 

650-850 

5B 

5.08 

2.0 

13.34 

5.25 

600-1200 

6B 

15.24 

6.0 

8.89 

3.5 

1000-1200 j 

7A 

15.24 

6.0 

8.89 

3.5 

800-1600 

7B 

7.62 

3.0 

8.89 

3.5 

800-1600 

8A 

2.54 

1.0 

14.61 

5.75 

1250-1500 

8B 

2.54 

1.0 

14.61 

5.75 

800-1600 

9B 

5.08 

2.0 

13.34 

5.25 

200-2000 


a Excitation signal decreases at a rate of 18 dB per octave outside indicated frequency range. 





TABLE V.-SUMMARY OF MEASURED DAMPING FACTORS 
USING POWER SPECTRAL DENSITY METHOD 


Run 

Mode order, m,n 

Frequency, f mn , Hz 

Damping, n 

□ a 

(a) High damping (spectral peak-to-valley ratio > 10 dB) 

2B 

2,1 

356 

0.0165 

1.08 

9B 

2,1 

359 

0.0141 

0.92 

3A 

3,1 

535 

0.0129 

1.12 

4B 

3,1 

538 

0.0128 

1.11 

9B 

3,1 

536 

0.0126 

1.10 

3B 

3,1 

539 

0.0144 

1.25 

5B 

4,1 

772 

0.0155 

0.89 

4B 

4,1 

772 

0.0196 

1.12 

7B 

2,3 

1176 

0.0133 

1.05 

7A 

3,3 

1346 

0.0107 

1.11 

9B 

3,3 

1348 

0.0105 

1.09 

(b) High damping (spectral peak-to-valley ratio < 10 dB) 

3B 

4,1 

780 

0.0216 

1.23 

5B 

3,2 

869 

0.0214 

1.89 

9B 

3,2 

865 

0.0120 

1.06 

7A 

2,3 

1185 

0.0155 

1.22 

7B 

3,3 

1340 

0.0126 

1.30 

(c) Low damping 

8B 

2,1 

364 

0.00098 

0.96 b 

4B 

3,1 

546 

0.00188 

1.10 

8B 

3,1 

546 

0.00193 

1.13 

8B 

1,2 

605 

0.00273 

1.35 c 

4B 

4,1 

789 

0.00174 

0.92 

8B 

4,1 

790 

0.00246 

1.29 

4B 

3,2 

889 

0.00121 

0.95 

8B 

3,2 

890 

0.00104 

0.82 

7A 

5,1 

1090 

0.00168 

1.32 

7B 

5,1 

1091 

0.00149 

1.17 

7A 

1,3 

1120 

0.00171 

0.91 

8B 

1,3 

1120 

0.00198 

1.06 

7A 

2,3 

1210 

0.00308 

0.98 

8B 

2,3 

1232 

0.00264 

0.84 

7A 

3,3 

1400 

0.00214 

1.10 

8B 

3,3 

1400 

0.00189 

0.98 


a Ratio of measured damping factor to datum damping factor. 

^Filter bandwith greater than one half of the resonance peak bandwidth. 
c Spectral peak-to-valley ratio < 10 dB. 




TABLE VI.— SUMMARY OF MEASURED DAMPING FACTORS USING 
SINGLE-SIDED FOURIER TRANSFORM METHOD 


Run 

Mode order, m,n 

Frequency, f m n , Hz 

Damping, 

B r a 

(a) High damping (spectral peak-to-valley ratio, > 10 dB) 

2B 

2,1 

357 

0.0156 

1.02 

9B 

2,1 

358 

0.0126 

0.82 

3B 

3,1 

537 

0.0139 

1.21 

4B 

3,1 

534 

0.0120 

1.04 

8B 

3,1 

536 

0.0123 

1.07 

9B 

3,1 

537 

0.0100 

0.87 

4B 

4,1 

775 

0.0188 

1.07 

4B 

4,1 

775 

0.0171 

0.98 

4B 

4,1 

770 

0.0152 

0.87 

8B 

4,1 

776 

0.0164 

0.94 

9B 

4,1 

773 

0.0150 

0.86 

7B 

2,3 

1179 

0.0126 

0.99 

8B 

2,3 

1194 

0.0148 

1.17 

7A 

3,3 

1347 

0.0105 

1.08 

8B 

3,3 

1354 

0.0080 

0.82 

(b) High damping (spectral peak-to-valley ratio < 10 dB) 

8B 

3,2 

862 

0.0183 

1.62 

7A 

5,1 

970 

0.0354 

1.77 

8B 

5,1 

970 

0.0178 

0.89 

7A 

4,2 

1056 

0.0205 

1.26 

7A 

1,3 

1091 

0.0214 

1.69 

8B 

1,3 

1085 

0.0120 

0.94 

7B 

3,3 

1340 

0.0088 

0.91 

(c) Low damping 

8B 

2,1 

364 

0.00174 

1.70 b 

8B 

3,1 

546 

0.00227 

1.34 b 

8B 

4,1 

790 

0.00216 

1.14 

8B 

3,2 

890 

0.00134 

1.05 

8B 

5,1 

1099 

0.00213 

1.68 

8B 

1,3 

1117 

0.00279 

1.48 

8B 

3,3 

1398 

0.00185 

0.96 


a Ratio of measured damping factor to datum damping factor. 
b F ilter bandwidth too wide (Afp > 1/4 Af^) 
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g'^= damping factor estimated when 
background noise is present 
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Figure 3.— Effect of Off-Resonant Contribution on Measurement of Damping Factor 
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Displacement power spectral density ratio (resonance peak to background noise valley) 



Displacement power spectral density ratio (resonance peak to background noise valley) 
Figure 4. —Correction to Damping Factor (Background Noise Effect) 
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Figure 8.— Effect of Lag Window on Parameter — 
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Computer software 
for SSFT 


Figure 10.— Flow Chart of Digital Data Reduction 












96 


FIGURE 12.— WORK ELEMENTS AND FLOW CHART FOR DEMONSTRATION EXPERIMENTS 

















Mode (1,1) 



Mode (2,1) 



Mode (3.2) Mode (5,1) 

FIGURE 14.-NODAL LINE PATTERNS FOR TEST PANEL 
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Mode (2,3) 



Mode (3,3) 



Mode (4,3) 



Mode (6,2) 



Mode (7,1) Mode (5,3) 


FIGURE 14.-CONCLUDED 
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Response phase angle 

(with respect to excitation), 0 Response amplitude 






(f) Vector Diagram of SSFT 


(c) Kennedy-Pancu Vector Diagram 

Figure 15. — Three Schemes for Calculating Damping Factor 
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1184.0 Hz 



Figure 17. — Vector Diagram for Mode (2,3), High Damping 
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Figure 18.— Vector Diagram for Mode (2,1), High Damping 
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Displacement power spectral density (arbitrary level) 
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Displacement power spectral density (arbitrary scale) 
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Displacement power spectral density (arbitrary scale) 



'9 















Figure 27. — Vector Diagram of F m ( to) for Modes (5, 1), (4,2), and ( 1,3), Run 7 A, High Damping 
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Frequency, Hz 

Figure 28.— Change in Arc Length With Frequency 


